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In  this  investigation  neutron  spectra  were 
measured  both  integrally  and  differentially  and  then 
compared  with  theoretical  calculations.   Measurements 
were  made  in  highly  absorbing  media  in  the  University 
of  Florida  Training  Reactor  (UFTR)  thermal  column,  in 
the  UFTR  core  and  in  the  subcritical  assembly  sitting  on 
top  of  the  graphite  pedestal  of  the  UFTR.   The  subcritical 
assembly  contained  some  typical  D20-moderated  natural 
uranium  lattices. 

Precise  knowledge  of  the  neutron  spectrum  in  a 
heterogeneous  medium  makes  it  possible  to  optimize  the 
neutron  economy.   This  is  true  because  fuel  elements  in 
a  reactor  have  varying  cross  sections  for  different  energy 
neutrons.   By  proper  choice  of  the  fuel-to-moderator  ratio 


ix 


and  the  geometry  and  composition  of  fuel  elements,  optimum 
use  of  neutrons  can  be  attained  for  maximum  fuel  burn-up 
with  the  highest  conversion  ratio. 

In  the  thermal  column  of  the  UFTR,  thermal  neutron 
spectra  were  measured  integrally  using  activation  detectors 
such  as  Au197,  Dy164,  Lu175,  Lu176  and  Pu239.   The  detectors 

Lu    and  Pu    have  low- lying  absorption  resonances, 

1P7        175 

whereas  Au    and  Lu    have  a  1/v  absorption  cross  section 

for  neutrons  below  1.0  ev.   The  detector  Dy    has  an 
absorption  cross  section  that  is  steeper  than  1/v.   The 
measurements  were  made  in  a  stainless  steel  rod,  a  natural 
uranium  slab  and  in  an  Al  can  filled  with  borated  water 
and  results  were  interpreted  using  Westcott's  formulation. 

In  a  natural  uranium  slab  the  effective  neutron 
temperature  change  from  the  center  to  the  outside  boundary 
of  the  slab  was  obtained  by  measuring  the  activity  ratios 
of  Lu176  and  Lu1   .   In  a  stainless  steel  rod,  the  spacial 
effective  neutron  temperature  change  was  measured  using 
Pu239  and  Lu176  activation  detectors  and  then  the  ratios 
of  activities  to  the  activities  of  1/v  absorbers  were 
obtained. 

The  same  type  of  measurements  were  made  in  an  Al 
can  filled  with  different  concentrations  B203#   These 
results  were  compared  with  THERM0S  calculations;  however, 
due  to  the  anisotropy  of  the  experimental  source  conditions, 
agreement  was  poor. 


The  effective  neutron  temperature  and  the  epithermal 
index,  Wt/T0|  were  measured  by  using  activation  detectors 
in  the  UFTR  core  between  the  two  fuel  boxes  in  the  North- 
South  center  line.   The  result  of  the  effective  neutron 
temperature  measurement  was  compared  with  the  open-beam 
differential  spectrum  and  found  to  be  in  good  agreement. 

A  crystal  diffraction  spectrometer  for  differential 
spectrum  measurements  in  the  subcritical  assembly  was  used. 
For  these  measurements,  the  total  efficiency  of  the  machine 
for  different  Bragg' 3  angle  was  first  obtained  by  taking  a 
neutron  beam  with  known  spectrum  out  of  the  subcritical 
tank.   The  tank,  24  in.  in  diameter  and  60  in.  long,  was 
filled  with  D20  and  rested  on  top  of  the  graphite  pedestal 
of  the  UFTR.   The  neutron  spectrum  out  of  this  medium  was 
measured  at  two  different  temperatures.   The  ratios  of  the 
two  spectra  were  plotted  with  the  ratios  of  the  calculations 
and  found  to  be  in  good  agreement. 

The  neutron  spectrum  out  of  the  subcritical  assembly 
was  measured  for  two  different  lattice  pitches  using  Mark  I 
and  Mark  V-B  natural  U  fuel  elements  in  D20.   A  neutron  beam 
was  extracted  from  the  cell  boundary  and  from  the  center  of 
the  cell  for  each  of  the  three  lattice  arrangements.   The 
results  were  compared  with  the  theoretical  spectra  calculated 
with  the  THERM0S  code  and  found  to  be  in  excellent  agreement. 

Integral  measurements  were  also  made  in  the  sub- 
critical  using  Lu176  and  Au197  detectors.   The  experimental 


activation  measurements  were  compared  with  the  THERM0S 
calculations  and  were  also  found  to  be  in  excellent 
agreement. 

Finally,  the  effective  neutron  temperature  change 
in  a  unit  cell  for  four  U-D20  lattice  arrangements  were 
measured  differentially  and  integrally  and  compared  with 
the  theoretical  calculations.   The  results  were  in  good 
agreement. 


CHAPTER  I 

DESCRIPTION  OF  THE  PROBLEM 

Neutron  spectrum  measurements  in  heterogeneous 
media  are  important  for  several  reasons:   (1)  for 
evaluation  of  the  thermal  utilization  factor  f  and 
the  thermal  reproduction  factor  n;   (2)  for  the  cal- 
culation of  reactor  fuel  cycles;   (3)  for  comparison 
with  multigroup  calculations  of  the  neutron  spectrum; 
and  (4)  for  the  prediction  of  dose  rates  in  a  sample, 
which  depends  on  the  neutron  energy  distribution. 

Precise  knowledge  of  the  neutron  spectrum  in 
a  heterogeneous  medium  makes  it  possible  to  optimize 
the  neutron  economy.   This  is  true  because  fuel 
elements  in  a  reactor  have  varying  cross  sections  for 
different  energy  neutrons.   By  proper  choice  of  the 
fuel-to-moderator  ratio  and  the  geometry  and  composi- 
tion of  fuel  elements,  optimum  use  of  neutrons  can 
be  attained  for  maximum  fuel  burn-up  with  the  highest 
conversion  ratio.   In  many  cases,  by  controlling  the 
leakage  spectrum  from  a  reactor  core  to  a  desired  form, 
the  maximum  conversion  ratio  in  the  blanket  material 


will  be  realized.  This  follows  because  U238  and  Th232 
which  are  the  usual  fertile  blanket  material  for  nuclear 
reactors,  have  high  resonances  for  neutron  capture  at 
6.8  ev  and  22-24  ev,  respectively.   In  order  to  pre- 
dict the  spacial  and  spectral  distribution  of  neutrons 
theoretically,  the  behavior  of  neutrons  in  the  hetero- 
geneous media  must  be  well  understood. 

In  this  investigation,  the  neutron  spectrum  was 
measured  both  integrally  (1)  and  differentially  ( 2_) ,  and 
then  compared  with  theoretical  calculations  (3,  4,  5). 

Measurements  were  carried  out  in  a  highly 
absorbing  medium  in  the  University  of  Florida  Training 
Reactor  (UFTR)  thermal  column,  in  the  UFTR  core,  and  in 
the  subcritical  assembly  sitting  on  top  of  the  graphite 
pedestal  of  the  UFTR.   The  subcritical  assembly  con- 
tained some  typical  D20-moderated  natural  uranium 
lattices. 


Underlined  numbers  in  parentheses  refer  to  the  List 
of  References. 


CHAPTER  II 

MEASUREMENT  OF  NEUTRON  SPECTRA 
Introduction 

As  early  as  19  36,  neutron  diffraction  by 
crystalline  material  was  suggested  by  Elsasser  (6_) 
and  experimentally  demonstrated  by  Halban  and  Preiswerk  (.7) 
and  by  Mitchell  and  Powers  (£) .   Using  Ra-Be  sources, 
these  experiments  proved  that  neutrons  were  diffracted, 
but  they  were  by  no  means  able  to  obtain  monochromatic 
neutrons  due  to  the  weakness  of  the  sources. 

Soon  after  the  development  of  the  nuclear  reactor, 
neutrons  were  available  in  great  quantity,  and  neutron 
diffraction  spectrometers  were  built  at  Argonne  (9_,  ICO  , 
Oak  Ridge  (11),  Chalk  River,  Canada  ( 12)  ,  Harwell, 
England  (13_)  and  elsewhere.   Also,  at  the  same  time,  time- 
of-flight  techniques  were  developed  (14,  15).   About  the 
same  time  integral  methods  of  neutron  spectrum  analysis 
were  described  by  Westcott  ( 16_)  ,  Campbell,  et  al.  (17) 
and  Bigham  and  Tunnicliffe  (18).   Trice  (19)  and  Nisle  (^0) 
also  shed  light  on  this  subject. 


An  early  theoretical  treatment  of  the  spectral 
hardening  was  carried  out  by  Plass  ( 2JL) ,  who  assumed 
neutrons  of  each  energy  diffused  independently  without 
energy  interchange.   Later  models  were  developed  to 
calculate  neutron  spectra  by  the  use  of  free  hydrogen 
and  bound  hydrogen  scattering  kernels  (22,  23). 

One  of  the  better  methods  for  calculating  the 
space  and  energy  dependent  flux  in  a  unit  cell  was 
developed  by  Honeck  (24,  2£)  for  H20  and  D20-moderated 
systems;  it  is  described  in  Appendix  II.   Recently 
at  KAPL,  Goldman  (26)    constructed  a  scattering  law 
based  on  Nelkin's  model,  and  agreement  with  the 
experimental  results  using  (CH^  ia  excellent.   At 
General  Atomics,  Young  C 2_7_)  has  measured  neutron  spectra 
using  a  pulsed  high-current  electron  linear  accelerator 
in  H20,  CH2 ,  and  CgHg  moderators  poisoned  with  B,  Sm,  Er 
Gd  and  also  D20  ( 2_8) .   Measured  spectra  were  compared 
with  DSN  (5)    transport-theory  calculations  utilizing 
the  bound  hydrogen  scattering  model  for  water.   Earlier, 
Beyster  (2£)  measured  neutron  spectra  in  pure  and 
poisoned  H20,  <CH2>n»  and  ZrH  using  the  pulsing  technique 
and  compared  it  with  DSN  calculations. 


Integral  Method 

The  measurement  of  neutron  energy  spectra  using 
the  time-of-flight  technique  or  a  crystal  spectrometer 
requires  the  extraction  of  a  neutron  beam  from  a 
reactor  lattice,  which  is  difficult  to  do  for  more 
than  a  few  representative  positions.   However,  by 
using  activation  detectors,  spectrum  changes  throughout 
the  cell  may  be  mapped  out  (i.e.,  good  spacial 
resolution  is  obtained  with  minimal  cell  perturbation). 

The  usual  experimental  technique  for  observing 
the  spectral  change  in  a  unit  cell  of  a  reactor  involves 
mapping  the  activity  ratio  of  resonance  and  1/v  cross 
section  detectors.   For  this  type  of  measurement,  various 
types  of  detectors  are  summarized  in  Table  2.1.   In 
Figure  2.1  the  cross  sections  of  the  1/v  absorber  Lu1  t 
and  the  resonance  absorber  Lu176(30)  are  shown.   Also 
included  in  the  figure  are  the  cross  sections  for  Dy 
( 31)  and  Pu2    (32).   A  Maxwellian  neutron  distribution 
in  the  thermal  energy  region  with  an  average  temperature 
of  29  3°K  is  also  shown  in  the  same  figure. 

By  inspection  of  Figure  2.1,  it  is  seen  that, 
as  the  Maxwellian  distribution  shifts  to  higher  energies, 
the  ratio  of  the  activities  of  Lu176  and  Lu175  will 
increase  since  the  reaction  rate  is  proportional  to  the 
integral  of  the  neutron  flux  times  the  absorption  cross 
section  ( 3 3) l 


Figure  2.1.   Activation  Cross  Sections 
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In  well-moderated  media  with  low  absorption, 
(e.g.,  D20,  H20,  C,  Be)  far  away  from  neutron  sources 
and  boundaries,  the  spectrum  approaches  a  Maxwellian 
distribution  with  the  average  temperature  corresponding 
to  the  temperature  of  the  medium.   When  these  well- 
thermalized  neutrons  enter  a  highly  absorbing  medium 
(e.g.,  a  fuel  element  which  has  a  1/v  absorption  cross 
section),  lower  energy  neutrons  are  preferentially 
absorbed,  resulting  in  a  "hardening"  of  the  energy 
spectrum  within  the  fuel. 

The  activation  results  were  interpreted  by 
using  the  Westcott  formulation  (34.) .   This  method  applies 
to  well-moderated  systems  such  as  the  UFTR.   In  a  reactor 
spectrum,  effective  cross  sections  are  given  as: 


wh< 


o 


8      =   ao  (g+rs)  (2.1) 


effective  cross  section 

cross  section  at  v  =  2200  m/sec 


r   ■   relative  intensity  of  the  slowing 

down  spectrum  (or  epithermal  index) 

g,s   =   functions  of  the  neutron  temperature,  T 


The  values  of  g  and  s  are  tabulated  ( 35_)  for  various 
nuclides.   For  1/v  absorbers  g  =  1  and.s  =  0.   Using 
this  notation,  the  Cd  ratio  for  a  very  thin  detector 
is  given  as: 


'cd 


g  +  rs 


.  gr  /  T 
rs  +  — A/  — 


(2.2) 


where  1/K  is  approximately  the  density  fraction  of 
epithermal  neutrons  transmitted  by  the  Cd  (,3k).      The 
factor  1/K  varies  with  the  thickness  of  the  Cd  filter 
and  is  tabulated  in  the  literature  ( 3»Q .   For  1/v 
absorbers  g  =  1  and  equation  (2.2)  becomes 


1  ♦  rs 


<cd 


■£ 


K  v 


(2.3) 


Rearranging  this  equation  gives: 


*V: 


s(T  /T)    (R„h-D  ♦  Rrd/K 


cd 


(2.i+) 


Therefore,  by  using  very  thin  detectors,  the  factor 


r  v —  can  be  measured, 
To 


For  detectors  with  appreciable  self  absorption, 
equation  (2.4)  becomes  (18) 
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'th 


(1  -  R 


cd 


(FR   -1) 
cd 


Gr  8q 


R  .  (1/K  -W) 
cd 


(2.5) 


where 


!th   =   thermal  neutron  self  shielding  factor, 
G   =   resonance  self  shielding  factor, 


F  ■   Cd  transmission  for  resonance  flux, 

W   =   fraction  of  resonance  activation  below 
Cd  cutoff, 

h   a   thermal  transmission  of  Cd  filter. 


In  order  to  measure  the  effective  temperature 
of  the  thermal  neutrons,  the  activity  ratios  of  Lu 
and  Lu176m  were  measured.   The  activity  ratios  of  these 
isotopes  at  position  x  with  respect  to  a  standard 
position  s  is  expressed  as 


A177 

X 

(g  +  rs)176 

X 

A176m 

_(g  +  rs)175_ 

YX 

Al77  " 

A176m 

s 

(g  +  rs)176 

Ys 

_(g  +  rs)175^ 

(2.6) 


*   Equation  (2.6)  is  for  very  dilute  detectors,  where 
resonance  and  thermal  self  absorption  is  negligible, 
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By  assuming  Gth  =  1  and  G  s  1,  and  solving  for  Y  , 
the  following  relationship  is  obtained: 


(g  +  rs  Gr) 


176 


Y   =  RY 
x     s 


(g  +  rs  G^) 


175 


(2.7) 


For  r  =  0  (i.e.,  no  epithermal  neutrons)  equation  (2.7) 
reduces  to 


176 


Yx  =  [g    ]x  =  RYS  =  R[gi76^ 


(2.8) 


The  effective  neutron  temperature  for  K176  (temperature 
index  of  Lu176)  is  obtained  from  Figure  2.2  (36). 

If  r  =  0  the  [(g  +  rs)176]^  must  be  evaluated. 
For  Lu175,  g  r  1  and  rVT/T0  is  a  measurable  value  with 
dilute  foils. 


sVT0/T  Gr I 
Jx 


a„  G, 


175 


'*  This  equation  takes  into  account  the  self-shielding 
of  resonance  neutrons. 


12 


Mrt  HW-«4««« 


4.0 

/*** 

1.3 

/ 

/  - 

J.O 

/ 

2.3 

/ 

J                       O    EXPERIMENTAL    VN.UU 

J                         •    CALCULATED  VALUE! 

a.o 

/ 

/                                  COUNTZRO    •TATIiTICf   1    1/M 

1.5 

•          / 

1.0 

200   300   400   500    600   700    900 


Figure  2.2.   Calibration  Curve  for.  g176 
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is  also  determined  experimentally.   Then  using  the  values 
°^  ^176  anc*  si+  given  by  Westcott  ( 3_5^  ,  corrected  values 
°^  ^176  and  s4  are  obtained  by  trial  and  error.   The 
effective  neutron  temperature  is  found  from  the  value 
of  8176  (see  Figure  2.2). 

An  alternative  method  of  evaluating  r"VT/T   is 
to  irradiate  two  detectors  and  express  their  activation 
ratios  as  follows : 


gGth  ♦  rs  Gr)' 
_(gGth  +  rs  Gr)] 


UA 


E 


(gGth  ♦  rs  G  )a  1 


(gGt,  ♦  rs  G  )] 


(2.9) 


or  rearranging 


CgGth  +  WtTt^  s^  Gj' 


o  o  r  x 


X  a  R'  X 

X  £ 


Cg61.h  ♦  rV  T/T 


Jth 


g  r 

o  o  r  x 


(2.10) 


Xs  is  determined  from  the  Cd  ratio  of  one  of  the  detectors, 
Solving  for  r"V  T/TQ  in  equation  (2.10),  the  following  is 
obtained: 
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R'    Xs(gGth)b    -    (gGth)' 


<rVT/TQ)x 


(s    G    )a  -   R'X   (s    G    )b 

or  s      o  r  (2.11) 


Differential  Method 

Time-of-Flight 

The  time-of-flight  method  utilizes  a  mechanical 
rotor  with  narrow  slits  called  a  chopper  (37)  or  a 
pulsed  accelerator  to  produce  bursts  of  neutrons  with 
a  time  duration  depending  on  the  speed  of  the  rotor  or 
accelerator  pulsing  system.   The  bursts  of  neutrons  travel 
through  an  evacuated  tube  to  a  bank  of  detectors .   The 
energies  of  neutrons  are  determined  by  electronically 
measuring  the  time  it  takes  for  neutrons  to  travel  from 
the  source  to  the  detectors. 

Usually  choppers  are  classified  as  slow  and 
fast.   Slow  choppers  give  bursts  of  about  20  w  sec. 
duration  whereas  the  burst  of  a  fast  chopper  is  on 
the  order  of  1  p  sec.   Unlike  crystal  spectrometers 
which  produce  monoenergetic  neutrons,  choppers  select 
certain  velocity  neutrons  out  of  heterogeneous  beams. 
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The  energy  resolution  of  a  chopper  is  described 
by  the  equation  ( 38) 


&      =  -0.0279^  E  ♦  2±» 
E  mm  (2.12) 


where 

AE  =  uncertainty  in  energy, 

E  =  energy  of  the  neutron, 

At  =  uncertainty  in  flight  time  per  meter, 

m  =  flight  path,  meters, 

Am  =  uncertainty  in  flight  path. 

For  a  fast  chopper,  the  second  term  is  negligible 
and  the  resolution  reduces  to  the  first  term.   Better 
energy  resolution  is  obtained  with  longer  flight  paths. 
A  list  of  resolutions  together  with  flight  paths  for 
different  choppers  is  listed  by  Anderson  ( 38) .   Multi- 
channel analysis   determines  the  intensities  of  neutrons 
in  different  energy  intervals. 

Slow  choppers  can  be  constructed  much  more 
easily  than  the  fast  choppers.   This  is  due  to  the  fact 
that  for  the  former,  chopping  is  accomplished  with  thin 
layers  of  Cd,  whereas  it  takes  many  inches  of  plastic  and 
steel  for  fast  choppers  (  39 ) . 
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An  alternate  method  for  measuring  neutron  energy 
spectra  (40)  is  the  time-of-flight  method  with  a  pulsed 
source.   The  advantages  of  this  method  are:   (1)  the 
fuel  elements  get  low  irradiation,   (2)  the  rotor  cut- 
off function  does  not  have  to  be  determined,  and   (3) 
by  the  use  of  the  chopper  and  pulsed  source,  the  neutron 
spectrum  as  a  function  of  slowing  down  time  can  be 
determined. 

Crystal  Spectrometer 

The  use  of  crystal  diffraction  spectrometers  for 
detecting  monochromatic  neutrons  has  been  investigated 
by  many  authors  (9^,  1£,  ¥1-47).   Diffraction  is  a 
scattering  process.   A  regularly  arranged  series  of  atoms 
of  a  crystal  will  scatter  the  neutron  waves  in  all 
directions,  but  only  in  particular  directions  will 
scattered  waves  be  in  phase  and  reinforce  each  other  to 
form  a  diffracted  beam.   The  atoms  of  a  crystal  are 
arranged  periodically  in  parallel  planes,  so  that  in 
general  scattered  waves  are  out  of  phase  except  in  a 
few  directions  where  reinforcement  takes  place. 

The  diffraction  process  applies  to  X-rays  as 
well  as  to  neutrons.   The  fundamental  difference  between 
the  two  is  that  X-rays  are  scattered  by  the  orbiting 
electrons,  where  as  the  neutrons  are  scattered  by  the 
nuclei  of  the  atom.   For  this  reason,  the  X-ray 
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scattering  amplitude  is  proportional  to  the  atomic 
number  of  the  element,  while  the  neutron  scattering 
amplitude  shows  a  relatively  srilall  variation  with 
atomic  number. 

When  neutrons  from  a  collimator  impinge  on  a 
crystal  plane  of  a  single  crystal  at  an  angle,  9,  only 
neutrons  of  one  energy  are  diffracted  in  the  direction 
of  the  angle  26.   Therefore,  the  crystal  and  the 
detector  must  maintain  an  angle  ratio  of  one-half  in 
order  to  obtain  monoenergetic  neutrons. 

The  diffraction  of  the  neutron  beam  depends  on 
the  lattice  spacing  in  the  crystal.   This  space  varies 
depending  on  the  crystal  material  and  on  the  crystal 
axis  along  which  it  is  cut.   The  coherent  scattering  of 
neutrons  from  the  nuclei  of  atoms  in  a  single  crystal 
to  produce  a  monochromatic  neutron  beam  is  governed  by 
the  familiar  Bragg' s  relation: 

n  X   =   2d  sin  6  (2.13) 

where 

d  =  lattice  spacing, 

n  =  order, 

A  =  wave  length, 

9  =  glancing  angle. 

In  energy  terms,  this  relationship  becomes: 
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E(ev)  =  kn2/sin2  9 

where  k  is  a  constant  for  each  lattice  spacing. 
Table  2.2  gives  a  list  of  crystals  with  k  and  the 
detectable  energy  at  1° . 


(2.14) 


TABLE  2.2 


TABLE  OF  k  AND  ENERGY  FOR 
DIFFERENT  CRYSTALS  AT  1° 


Crystals 


k  x  10°  (ev) 


Enerj 

?y  C 

ev) 

15, 

,40 

8, 

.42 

12, 

,28 

6, 

.27 

12, 

,25 

Cu  (111) 

NaCl  (200) 

Al  (111) 

Ge  (111) 

LiF  (111) 


4.71 
2.58 
3.76 
1.92 
3.81 


From  equation  (2.14)  it  can  be  seen  that  the 

higher  order  interference  will  be  due  to  neutrons  with 

2 
energy  n'  times  the  lowest  order.   Also,  there  is  inter- 
ference due  to  the  mosaic  spread  of  the  crystal. 
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Mosaic  Structure 

In  order  to  know  the  energy  of  diffracted 
neutrons  at  different  angles,  it  is  necessary  to  know 
the  distance  between  the  atomic  planes  of  the  crystal. 
The  distance  depends  on  the  direction  in  which  the 
crystal  is  cut  and  also  on  the  crystal  material.   If 
the  crystal  were  perfect,  neutrons  of  a  single  energy 
would  be  diffracted  at  any  given  angle.   However,  all 
real  single  crystals  have  a  mosaic  structure,  that  is, 
they  possess  structural  imperfections  which  break  up 
the  arrangement  of  atoms  into  a  number  of  small  blocks , 
each  slightly  disoriented  from  one  another.   This 
mosaic  structure  causes  the  diffracted  neutrons  to 
have  a  spread  of  energies  at  any  angle. 

Higher  Order  Diffraction 

Also,  according  to  equation  (2.14),  there  are 
higher  order  neutrons  diffracted  at  each  angle.   The 
ratio  of  the  second  order  to  first  order  neutrons 
detected  in  the  beam  can  be  written  as  (^jO  : 

f(E)  =  4(E)  e(E)  «5(E)ns2  w(E)n=2  p(E)n=2 
where 


(2.15) 


6(E) 
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$(E)  =  lilEl  U(E)  is  assumed  to  be  1/E 

*(E)  for  energy  greater  than  0.4  ev) 

,n    e(4E) 

c^l;  -  (ratio  of  the  detector  efficiency 

e(E)  for  second  and  first  order  neutrons) 


AE(n=2,  HE)     1  (4E)3/2  cos[9(4E,  n=2)]A8(4E) 
n=2    AE(n=l,  E)      2   (E)3/2  cosC6(E,  n=l)]  A6(E) 


where 


(2.16) 


w(HE),n=2) 

w(E)n  =  2   =  sl    (w(E)  is  the  fraction 

w(E,  n=l)        of  reflected  beam  trans- 
mitted by  collimating 
system) 


r(4E)n_? 
p(E)  _0   =   2-i         (r(E)  is  reflectivity 


n=2 


of  the  crystal  plane 
used) 


In  equation  (2.16)   cos[6(4E,  n=2)]  =  cos  [8(E,  n=l)]. 

and  let  o  =  f  which  reduces  the  energy  band 

A6(E) 
transmitted  to  the  detector  to 


°<E)   o   =  a    a  (2.17) 

n=z      o 
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The  reflectivity  of  a  crystal  for  neutrons  is 
given  by  equation  (11)  of  Holm  (^9.)  for  the  Laue  Case 


_A   -    tan2e  (-lP*1 
r(E)  =  exp(-2M)  v/2    ne     £   Y  


where 


W       j  3 


exp(-2M)   ■   Debye-Waller  temperature 
correction  factor 

n   =   "mosaic  spread"  of  the  crystal 

A  =   product  of  the  linear  abosrption 
coefficient  and  the  path  length 
in  crystal 


(0     >in    n3 


(2.18) 


t   =   thickness  of  crystal 

N   =   reciprocal  of  the  unit  cell 

F  =   crystal  structure  factor 

The  Debye-Waller  temperature  correction  factor 
accounts  for  the  reduction  in  intensity  of  the  Bragg 
reflection  due  to  the  thermal  motion  of  the  atoms  in 
the  lattice,  where 

M  =  ^T^  (~)2  [1/4  +  (T/®)2  Q((H)/T)]        (2.19) 
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with 

.z   xdx 


Q(Z)   =   / 


o        ex_i  (2.20) 

Here, 

(h)   =   Debye  or  characteristic  temperature 
of  the  crystal 

h,k   s   Planck's  and  Boltzmann's  constants 

m  *   nuclear  mass 

T   s   temperature  of  the  crystal  in  degrees 
Kelvin 

Since  the  factors  in  equation  (2.15)  cannot  be 
calculated  accurately,  it  is  easier  to  obtain  an  open 
beam  count  rate  of  E  and  HE.   The  ratio  is: 


count  (HE) 
f*(E)  =  — 


count 


—  .  ♦(E)e(E),(E)nrlM(E,n=ip(E)n=i 


(2.21) 


Equation  (2.21)  is  a  good  approximation  because  the 
second  order  contamination  is  usually  on  the  order  of 
a  few  per  cent,  unless  the  first  order  energy  is  below 
the  Maxwellian  peak.   Now,  o(E)nsl  is 


AE(n*l,  HE)    HE3/2cos[6(HE,  nsl)]A9(HE) 
o(E)nsl  =  . 


AE(n=l,  E)       E  cos[9(E,  n»l)]Ae(E)         (2.22) 
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where 

cosCeC+E,  n=l)] 

6  .  

cos[e(E,  n*l)] 


o(E)n=1   =   2ooaS  (2.23) 


Now  the  f(E)  term  in  terms  of  f*(E)  is  expressed  as: 
r(UE)n=2 

f(E)  r  f*(E)  


r(4E)nsl  2Bw(E)nsl  (2>24) 

Since   Bw(E)nsl  ■  j,f  then 


r(«'„=2 


f(E)  =  f*(E)  1/2  =  f*(E)  R(4E) 

r(4E)n=l 


(2.25) 


Therefore,  by  using  the  calculated  value  of  R(HE),  f(E) 
can  be  estimated. 

The  Renninger  Effect 

When  a  neutron  spectrum  was  measured  using  a 
crystal  spectrometer  and  corrections  were  made,  many 
large  dips  were  found  (2_.  50-53)  instead  of  the  smooth 
spectrum  characteristic  of  a  Maxwellian  distribution 
of  neutron  velocities.   These  fluctuations  are  due  to 
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Bragg  reflection  of  neutrons  from  reflecting  planes 

other  than  that  used  to  obtain  the  monochromatic  beam. 

This  reflection  is  best  explained  by  the  use  of  the 

reciprocal  lattice  and  sphere  reflection  concepts  (54) 

and  the  vector  notation  as  outlined  in  Appendix  III. 

Figure  2.3  demonstrates  multiple  reflection  by 

means  of  the  sphere  of  reflection,  where  the  radius  of 

the  sphere  is  1/X;  it  is  centered  on  the  crystal  at  C, 

and  passes  through  the  origin  of  the  reciprocal  lattice 

at  0.   For  an  incident  direction  CO,  reflection  occurs 

in  the  direction  CP  if  the  reciprocal  lattice  vector 

r' ,  terminates  at  P  on  the  sphere.   Suppose  a  second 

reciprocal  lattice  vector  r^  terminates  at  Q  on  the 

sphere,  then  reflection  of  the  same  wave  length  occurs 

in  the  direction  CQ.   Taking  Q  as  the  origin  of  the 

reciprocal  lattice  and  CQ  as  the  incident  direction, 

further  reflection  corresponding  to  the  vector  r'  i8 

possible.   The  double  reflection  by  way  of  r'  and  r' 

gives  a  resultant  in  the  same  direction  and  at  the  same 

wave  length  as  a  reflection  through  r' .   Since  r'  =  r' 

-► 
+  r^  the  Miller  indices  of  the  planes  are  related  by: 

hx  =  h2  +  h3  ,  k1   =  k2  +  k3,  i±   =  *2  +  a  (2.26) 

The  double  reflection  can  occur  when  any  reciprocal  lattice 
point  Q  is  on   the  sphere,  where  Q  is  not  necessarily 
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Reflected  Direction 


q        Incident  Direction 

< 


Figure  2.3.   Sphere  of  Reflection 
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located  coplanar  with  C,  0,  and  P.   Finally,  it  should 
be  noted  that  when  multiple  Bragg  reflection  occurs 
for  first  order  reflection,  it  also  occurs  for  higher 
order  reflections  at  the  same  angle  of  incidence. 

Reflectivities 

Since  reflectivities  for  NaCl  and  Be  have  been 
calculated  ( 4_9) ,  a  NaCl  crystal  cut  along  the  (200) 
plane  was  used.   The  disadvantage  of  this  crystal  is 
that  the  second  order  contamination  becomes  very  large 
at  low  energies  (e.g.,  =  50  per  cent  at  E  =  .025  ev) . 
In  order  to  reduce  this  second  order  contamination, 
a  crystal  with  suppressed  second  order  reflectivity 
must  be  used;  LiF  (111)  is  such  a  crystal.   The 
suppression  of  second  order  reflection  is  due  to  the 
opposite  sign  of  scattering  from  successive  planes  (12) . 

There  are  some  disadvantages  in  using  LiF 
crystals.   These  are:   the  lattice  parameters  are 
uncertain  (_55)  ,  the  diffuse  background  scattering  is 
fairly  high,  single  crystals  are  hard  to  obtain  and 
suppression  of  second  order  is  not  very  great  [F  (111)/ 
F  (222)  =  1.97]. 

Use  of  Ge  and  Si  crystals ,  which  have  diamond 
structure, is  very  favorable  because,  for  the  (111)  plane, 
the  second  order  (222)  vanishes  completely.   In  fact, 
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it  has  been  shown  (SO)    that  at  0.015  ev,  the  higher  order 
contamination  of  a  Ge  (111)  crystal  is  only  2.5  per  cent. 
Use  of  a  quartz  filter  also  reduces  the  higher  order 
contamination  considerably  (57). 

Transmission  of  Reflected  Beam  Through  Filters 

Transmission  of  a  neutron  beam  through  a  filter 
can  be  expressed  as: 


t  =  iA  ♦  V2  ♦  t3x3  ♦  ...  (2>27) 


where 


transmission  of  neutrons  at 
energy  E. 

fraction  of  neutrons  at  energy 


The  sum  of  the  fractions  of  neutrons  must  add  up  to 
one  (EX^  =1).   By  use  of  absorbers  such  as  Au,  Cd, 
Sm,  Ag,  Rh,  In  etc.,  neutrons  with  certain  energies 
can  be  suppressed.   These  absorbers  have  a  large  resonance 
and/or  a  1/v  cross  section.   For  instance,  for  a  1/v 
absorber,  suppression  of  the  first  order  by  99  per  cent 
suppresses  the  second  order  by  90  per  cent.   Therefore, 
the  transmission  method  is  a  very  useful  one  for 
estimating  higher  order  contamination.   The  energy 
resolution  of  a  crystal  spectrometer  neglecting  the 
mosaic  spread  of  the  crystal  is  expressed  as  C+4): 
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R  =  —  =  2  Ae  cot  e 

E  (2.28) 


The  actual  resolution  is  determined  by  plotting  the  rocking 
curve  for  different  angles.   The  rocking  curve  is  ob- 
tained by  counting  diffracted  neutrons  with  a  stationary 
detector  and  by  rotating  or  rocking  the  crystal  to 
different  angles.   The  maximum  of  the  distribution 
corresponds  to  the  position  where  the  detector  and  crystal 
plane  are  aligned.   The  width  of  the  peak  at  half- 
maximum  is  equal  to  AE,  the  energy  spread,  of  the  system. 

The  main  advantage  of  a  crystal  spectrometer 
compared  to  the  time-of-f light  method  in  the  low  energy 
range  is  that  the  area  of  the  beam  at  the  detector  is 
quite  small.   A  diffraction  spectrometer  can  cover  the 
energy  range  from  0.01  to  10  ev,  whereas,  it  is 
difficult  to  cover  this  energy  range  with  one  chopper. 
The  disadvantages  of  a  crystal  spectrometer  are:   (1) 
the  reflectivity  varies  inversely  as  the  neutron  energy; 
(2)  the  contribution  of  the  higher  order  diffraction 
increases  below  the  Maxwellian  peak,  and;   (3)  the 
spectrometer  detector  is  closer  to  the  reactor  than  the 
time-of-f light  detectors,  resulting  in  a  higher  background. 


CHAPTER  III 

DESCRIPTION  OF  APPARATUS 

Crystal  Spectrometer 

Collimator 

In  collimating  neutrons  through  a  narrow  channel, 
the  intensity  decreases  very  rapidly.   Therefore,  it 
is  necessary  to  place  several  collimating  channels  side 
by  side,  where  each  channel  has  the  required  small 
angular  divergence.   This  type  of  collimator  is  called 
a  Sollar  collimator  (58) . 

A  Sollar  collimator  was  constructed  from  eight 
stainless  steel  strips  0.010  in.  thick,  1.50  in.  wide, 
and  33  in.  long.   Stainless  steel  rods  one-eighth  of 
an  inch  square  were  used  for  separating  the  strips. 
This  collimator  was  then  inserted  into  a  Cd-lined  slot 
1  in.  x  1.50  in.  x  33  in.  in  a  3  in.  diameter  wooden 
cylinder;  the  collimator  was  held  in  the  slot  by 
friction.   The  collimator  was  then  dipped  into  a 
boiling  bath  of  paraffin  wax,  (CH2)n,  dissolved  in 
benzene,  (CgHg).   This  procedure  gave  a  thin  coat  of 
(CH2)n  on  the  surface  of  the  collimating  plates  in 
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order  to  eliminate  the  secondary  scattering  from  the 
walls.   The  angular  divergence  of  this  collimator  is 
13.0  minutes;   one  section  of  the  collimator  is  shown 
in  Figure  3.1. 

Neutrons  entering  the  collimator  propagate  in 
four  different  ways:   (1)  direct  transmission   (2) 
mirror  reflection   (3)  transmission  through  the 
collimator  walls,  and   (4)  rescattering  by  the  walls. 
These  effects  are  now  considered  in  turn. 

For  direct  transmission  with  an  isotropic 
source  on  one  end,  consider  a  narrow  slit  with  thickness 
T,  width  W,  and  length  L,  as  shown  in  Figure  3.1. 

The  number  of  neutrons  emerging  from  the 
collimator  per  unit  width  .is  found  from  the  relation  (59) 


9  1      T 

N(neutrons/cm  sec)  =  /  

o      2*    L,2Sec2    e 


x      <|>o   TdlL      tan    6)  (3.1) 


where 

d(L   tan    9)  =      line. source   element 

♦o 

~  Td(L  tan  9)   =   source  strength  per  unit 
solid  angle 
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Figure  3.1.   Single  Collimator  Slit 
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Figure  3.2.   Size  of  the  Effective  Source 


32 


distance  from  emitting 
element  d(L  tan  8)  to  the 
receptor 


angle  from  the  emitting 
element  to  the  receptor  in 
radians 


On  integrating  equation  (3.1),  the  following  result  is 
obtained: 

T2 
N(neutrons/cm  sec)  =  — - —  <•>   e 

2»  L   °  (3.2) 

Dividing  both  sides  of  equation  (3.2)  by  thickness  T, 
a  relationship  between  the  direct  inlet  and  outlet 
scaler  flux  is  obtained  as  follows: 

o  ^*  9 

*   (neutrons/cm^  sec)  =    ° 

2ff  L  (3.3) 

The  intensity  of  the  neutron  beam  emerging  from 
the  collimator  would  be  a  triangular  function  if  the 
collimator  were  perfect: 

I  <♦)   =   I0  Cl  "  ♦/o]  C3.«0 

where 

K*)   =   intensity  as  a  function  of  angle  $ 

I0   =   maximum  intensity  in  the  center  of 
the  collimator  axis 
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and  the  angle  $  varies  from  -a  to  +o.   Due  to 
imperfections  in  the  collimator  and  the  finite  trans- 
mission of  neutrons  through  the  wall,  this  triangular 
function  tends  to  be  closer  to  Gaussian  distribution 
(60).   The  Gaussian  approximation  of  equation  (3.4)  is 


I  <♦)   =   IQ  exp  -U/a')2  (3.5) 


where 

1/2 


a 


=   a/(ln2) 


In  the  actual  experimental  set-up,  the  disk 
source  is  located  some  distance  below  the  collimator 
inlet;  therefore,  an  approximate  relationship  must  be 
derived  between  the  scaler  flux  coming  out  of  the 
collimator  and  the  scaler  flux  of  the  source.   Again, 
assume  the  source  is  isotropic  and  consider  only  those 
neutrons  which  are  emitted  directly.   Figure  3.2  shows 
Tg  and  Tw ,  the  width  and  the  thickness  of  the  source 
subtended  by  the  collimator  outlet  which  is  approximately 
rectangular  in  shape. 

The  expressions  for  T  and  W  are  as  follows: 

Ts   =   T+JIJLj/lJt 

Wg   =   W+2  [l2/lJ  W  (3.6) 


Now  imagine  that  the  collimator  is  of  thickness  T 

s » 

width  Wg  and  length  (L^L^.  Uging  equation  (3.3), 
the  relationship  for  the  outlet  scaler  flux  for  thii 
collimator  is: 


T  ♦  ( 

o  s  ° 

♦'(neutrons/cm^  sec)   =  


2L'  (3.7) 

Since  the  actual  area  of  the  outlet  collimating  slit 
is  T-W  rather  than  T^,  an  approximate  expression  for 
the  outlet  scaler  flux  using  equation  (3.5)  is  obtained 
as  follows : 


vTs    J  ,W-c*/«v2]  d* 


♦  (neutrons /cm2  sec)  =  $'  ■■"  i. 


/  sexp[-(*/a»  )2]  d* 


(3.8) 


where 


.   1/2(T   +  T) 


'1 


♦  „  =   Arc  sin  tTs/LsJ 


35 

The  total  reflection  of  neutrons  by  the  wall 
of  the  collimator  can  distort  the  neutron  spectrum. 
Such  a  reflection  was  first  reported  by  Fermi  and  Zinn 
(61).      The  critical  angle  for  such  reflection  is  (62). 


,     ,1/2 
MNa/»)  (3.9) 


where 

X   =   neutron  wave  length 

N   =   number  of  atoms /cm3  in  the  scattering 
medium 

"a"   =   average  coherent  scattering  length 

Equation  (3.9)  shows  that  the  critical  angle  is  a 
function  of  neutron  energy.   In  the  relationship  for 
the  critical  angle,  "a"  can  be  either  positive  or 
negative,  depending  on  the  scattering  material.   For 
instance,  hydrogen  has  a  negative  scattering  length, 
while  carbon  has  a  positive  value.   Experiments  have 
been  performed  using  neutron  mirrors  with  various  hydro- 
carbons (&y>    and  it  was  found  that  6   becomes  zero  for 
a  ratio  of  H/C  -  1.748,  i.e.,  the  negative  scattering 
length  of  H  is  equal  to  the  positive  value  of  C.   There- 
fore, if  the  collimator  walls  have  a  value  of  H/C 
larger  than  1.748,  total  reflection  cannot  occur. 
Accordingly,  the  surface  of  the  collimator  walls  was 
coated  with  a  thin  coat  of  paraffin  which,  having  an 
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H/C  ratio  of  approximately  2,  should  eliminate  the 
total  reflection  of  neutrons.   This  type  of  collimator 
has  been  constructed  previously  (6JO  with  a  thin  film 
of  polyethylene. 

The  probability  of  direct  transmission  through 
the  wall  of  the  collimator  without  absorption  or 
scattering  is  expressed  as  e"   tz  (  Where  E   i3  the 
total  cross  section  of  the  wall  material  and  z  is  the 
path  length  of  the  neutrons  through  the  walls.   Even 
though  the  absorption  of  steel  is  not  high,  this 
probability  is  very  small  since  any  neutrons  which  the 
steel  absorbs  will  pass  through  it  at  very  small  angles 
Thus,  the  neutron  path  length  in  steel  is  relatively 
long,  and  transmission  through  the  wall  is  negligible 
compared  to  the  direct  transmission. 

Calculations  have  been  made  (2)  for  the  trans- 
mission of  neutrons  by  rescattering  in  the  walls  of  a 
similar  collimator.   For  this  calculation,  this  con- 
tribution from  rescattering  was  found  to  be  less  than 
one-half  per  cent  of  the  total  transmitted  intensity 
with  negligible  effect  upon  the  angular  distribution. 

The  description  above  states  that  the  neutrons 
emerging  from  the  collimator  come  predominantly  from 
direct  transmission.   A  calculation  was  made  for  the 
ratio  of  <J>/*0  for  the  experimental  set-up  and  was  found 
to  be  equal  to  7.14  x  10"5. 
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Crystals 

NaCl  (200)*  and  LiF  (111)*  crystals  were  used 
for  these  experiments.   Dimensions  of  both  crystals 
were  1/2  in.  x  2  in.  x  2  in. 

Goniometer** 


The  goniometer  consists  of  a  sample  table  with 
a  crystal  holder,  an  arm  on  which  the  second  collimator 
and  detector  shield  are  attached,  and  the  gearing 
system.   The  gearing  system  is  in  a  cast  aluminum  box, 
and  they  can  be  moved  manually  or  electrically.   Figure 
3.3  shows  the  crystal  spectrometer  in  operation  on  top 
of  the  University  of  Florida  Training  Reactor  (UFTR). 
The  crystal  holder  is  made  from  a  rectangular  aluminum 
box  with  both  ends  open  in  which  the  crystal  fits. 
This  box  is  attached  to  the  crystal  table  by  means  of 
a  stainless  steel  rod  3/4  in.  diameter. 

The  detector  shield  is  attached  to  an  arm  with 
a  collimator  and  is  shown  in  detail  in  Figure  3.4.   The 
movement  of  the  crystal  holder  and  the  detector  arm  can 


*  These  crystals  were  obtained  from  The  Harshaw 
Chemical  Company,  1915  E.  97th  Street,  Cleveland 
6,  Ohio. 

**   Picker  Diplanar  Diffractometer,  Model  No.  3527. 


Figure  3.3.   Diffraction  Spectrometer  in  Operation 
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be  performed  separately  or  at  an  angle  ratio  of  1:2. 
The  1:2  movement  is  carried  out  at  very  high  precision 
by  means  of  worm  gears  in  the  machine.   Two  odometers 
are  attached  to  the  box  where  the  angular  positions  of 
the  crystal  and  the  detector  can  be  read  to  l/100th  of 
a  degree. 

Neutron  Detection  System 

The  detector  used  for  these  measurements  is 
2  1/32  in.  diameter,  96  per  cent  B10  enriched  BF3 
counter  with  a  pressure  of  76  cm  Hg.   The  length  of  the 
tube  is  15  3/4  in.  long  with  12  3/4  in.  of  sensitive 
length;  it  has  a  ceramic  end  window.*   The  center  of 
the  plateau  for  this  detector  was  found  to  be  2  800  volts 
using  the  existing  cables.   The  length  of  the  plateau 
was  about  400  v.  with  change  of  slope  of  one  per  cent 
per  100  volts. 

The  efficiency  of  the  detector  for  different 
energy  neutrons  was  not  calculated  because  it  was 
determined  experimentally.   Theoretically,  efficiency 
can  be  calculated  from  the  following  relationship  (44): 

E  =  1  -  exp  (-NXCE-1/2)  (3.10) 


*  This  detector  was  supplied  by  Reuter-Stokes 

Electronic  Components,  Inc.,  Model  No.  RSN-108S. 
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where 

N   =   number  of  atoms  per  cm3  in  the 
counter 

x   =   length  of  the  counter 

C   =   numerical  constant  which  gives  the 
1/v  slope  of  the  BlO(n  a)  L^7 
cross  section 

E   =   energy. 

With  the  BF  counter,  a  Radiation  Counter 
Laboratory  (RCL)  decade  scaler,  preamplifier  and  linear 
amplifier  (model  numbers  2032,  20200  and  20100, 
respectively)  were  used.   In  addition,  Atomic  Instru- 
ment Company  regulated  high  voltage  supply,  Model  319, 
was  used. 

The  Subcritical  Assembly 

A  subcritical  tank  made  from  6061  Al ,  24  in. 
diameter  and  60  in.  long  and  wrapped  with  .030  in. 
thick  Cd,  was  placed  on  top  of  the  graphite  pedestal 
of  the  UFTR.   A  grid  made  from  the  same  type  of  Al  was 
placed  in  the  tank,  so  that  the  fuel  elements  could  be 
arranged  in  a  hexagonal  pattern  with  lattice  pitches 
of  14.7  and  22  cm. 

Figure  3.5  shows  the  location  of  the  tank  in 
the  UFTR,  and  Figure  3.6  shows  the  top  view  of  the 
lattice  arrangements  in  the  subcritical  assembly. 
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Figure 


3.6.   Top  View  of  the  Subcritical  Lattice 
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The  plane  source  of  the  neutrons  entered  the  subcritical 
assembly  from  the  bottom  face  of  the  tank.   There  is 
approximately  4  in.  of  moderator  between  the  grid  and 
the  bottom  of  the  tank. 

Two  types  of  natural  uranium  fuel  elements  were 
used  for  these  experiments:   Mark  V-B  and  Mark  I.   Mark 
V-B  is  an  annular  type  of  element  2.684  in.  o.d.  and 
1.970  in.  i.d.  and  8  9/32  in.  long.   Mark  I  is  a  solid 
fuel  element  one  in.  in  diameter  and  8  3/8  in.  long. 
There  were  six  fuel  elements  in  each  tube.   The  moderator- 
to-fuel  volumetric  ratio  Vm/Vf  and  pitches  are  listed 
in  Table  3.1. 

TABLE  3.1 

SUBCRITICAL^  LATTICE  CONFIGURATIONS 


Pitch  (cm) 

Mark 

V-B 

Rods 

Mark  I 

Rods 

Vvf 

22.0 

7 

0 

23.76 

22.0 

7 

7 

18.04 

14.7 

13 

0 

12.34 

14.7 

13 

13 

9.25 

Some  of  the  parameters  of  the  subcritical 
assembly  lattices  using  the  BUCKSHCT  Code  (6_5)  were 
calculated  and  are  listed  in  Appendix  IV.   A  proposal 
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for  performing  experiments  in  the  subcritical  assembly 
utilizing  the  UFTR  as  a  plane  source  is  outlined  in 
Appendix  E.   'The  critical  rod  position  of  the  reactor 
'was  checked  after  the  installation  of  the  subcritical 
assembly  as  outlined  in  the  proposal.   No  coupling 
between  the  UFTR  and  the  subcritical  assembly  was  found. 

Activation  detectors  were  used  for  the  integral 
measurements  and  a  crystal  diffraction  spectrometer  for 
differential  measurements.   For  comparisons  of  experi- 
mental results  with  the  theory,  THERM0S  (4_)  and  CEPTR 
(3)  codes  were  used. 

Two  five  kilowatt  heaters  were  inserted  into 
the  tank  so  that  the  D20  could  be  heated,    A  constant 
temperature  in  the  tank  was  maintained  by  the  use  of  a 
temperature  control  unit  which  utilized  four  thermo- 
couples.  Also,  a  stirrer  helps  to  keep  a  uniform 
temperature  throughout  the  tank.   A  side  view  of  this 
set-up  in  the  tank  is  shown  in  Figure  3.7. 

The  hydrogen  content  of  the  D20  was  determined  by 
use  of  the  NMR  (Nuclear  Magnetic  Resonance)  equipment  in 
the  University  of  Florida  Chemistry  Department*  and  found 
to  be  0.661.05  per  cent. 


*     Thanks  are  due  Dr.  Wallace  Brey  of  the  University  of 

Florida  Chemistry  Department  for  making  this  determination 
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Figure  3.7.   Side  View  of  the  Moderator  Tank 


CHAPTER  IV 

RESULTS  OF  INTEGRAL  SPECTRUM  MEASUREMENTS 

Introduction 

The  first  integral  measurements  were  made  in 
the  UFTR  core  and  thermal  column  where  there  jere  few  fast 
neutrons  (66_) .   The  spacial  hardening  of  the  thermal 
neutrons  in  highly  absorbing  media  was  observed. 
Then,  the  spectral  hardening  in  a  subcritical  exponential 
facility  which  uses  the  graphite  pedestal  of  the  UFTR 
as  a  neutron  source  was  measured  spacially  using 
differential  and  integral  methods  and  results  were 
compared  with  theoretical  calculations  (4). 

For  integral  measurements  of  the  neutron  spectrum, 
various  types  of  activation  detectors  were  used.   All  of 
the  detectors  fell  into  two  categories:   either  they 
had  a  1/v  absorption  cross  section  or  they  had  a  large 
low-lying  absorption  resonance. 
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Measurements  in  the  Thermal  Column* 

Stainless  Steel  Rod 

Integral  measurements  were  made  in  a  stainless 

steel  rod  one  in.  diameter  and  12.5  in.  long  and  in  a 

uranium  slab  2  in.  x  2  in.  x  1/2  in.   Radial  holes 

0.050  in.  diameter  and  1/2  in.  deep  were  drilled  in 

the  stainless  steel  pin  holder  6.5  in.  from  the  end. 

This  was  done  by  cutting  the  stainless  steel  rod  into 

two  pieces.   These  holes  held  wire  detectors.   For  the 

uranium  slab,  detectors  were  sandwiched  in  between  the  U 

foils  which  were  2  in.  x  2  in.  x  .005  in.   Measurements 

were  made  2  4  in.  from  the  east  end  of  the  thermal  column 

of  the  UFTR.   At  this  position,  the  ratio  of  thermal  to 

epithermal  flux  is  approximately  100  (66).   Therefore, 

the  contribution  of  the  epithermal  neutrons  is  negligible 

and  equations.  8)  is  used  directly  to  measure  the  change 

in  effective  neutron  temperature. 

For  these  measurements,  Lu175,  Lu176,  Dy16\ 
.19  7         2  39 
Au    ,  and  Pu    detectors  were  used.   The  detectors 

were  fabricated  in  the  form  of  wires  0.0  30  in.  diameter 

containing  10  per  cent  by  weight  of  the  activant  in  an 

Al  matrix,  except  for  the  Au  wires,  which  were  0.005  in. 

in  diameter.   The  first  measurements  were  made  in  the 


*  Preliminary  results  given  in  BNL-719,  Vol.  II. 
pages  336-58  (1962). 
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stainless  steel  rod  using  natural  Lu-Al  wires.*   a  few 
hours  after  irradiation,  Lu1  m  and  Lu177  activities 
were  counted  five  times  using  a  scintillation  counter.** 
About  four  days  later,  after  the  Lu     activity  had 

decayed  out,  five  more  sets  of  counts  were  taken  of  the 

177 
Lu    activities.   The  total  count  rate  for  each  time 

was  at  least  10,000  counts.   The  procedure  for  obtain- 
ing the  average  normalized  activities  of  the  two 
isotopes  and  their  ratios  is  outlined  in  Appendix  A. 

Pu-Al  wires  were  calibrated  by  measuring  their 
natural  gamma  activities.   These  wires  were  then 

irradiated  in  a  stainless  steel  rod  and  gammas  from 

19  7 
the  fission  fragments  were  counted  (67).   The  Au 

and  Dy    wires  were  also  irradiated  in  the  stainless 
steel  pin  holder.   The  corrected  activities  from  all  the 
detectors  are  listed  in  Table  4.1.   The  activities  of 
the  1/v  abosrbers,  the  Lu1   and  the  Pu239  fission  pro- 
ducts were  each  fitted  to  the  Bessell  function,  I  (<K  r) 
and  are  plotted  on  Figure  4.1.   The  ratios  of  fitted 
curves  f or  (Lu   /i/v  absorber)  and  (Pu239/l/v  absorber)*** 


*   Natural  Lu  contains  2.6  per  cent  Lu^76  and 
97.4  per  cent  Lu175. 

**   Baird-Atomic,  Inc.,  Gamma  Spectrometer,  Model  No.  810, 

***   Activities  of  Pu239  means  activities  of  the  fission 
fragments  of  Pu239. 


O  Au198  act. 
A  Lu176m  act 
*  Lul77  act# 
°  Pu239  act> 


/lo(1.045r) 


Io(0.980r) 


Io(0.880r) 


Io(0.980r) 

- 

-Io(1.0«*5r) 

Io(0.880r) 

Io(0.980r) 


Rod 
Boundary 


0.5 


1.5 


1.0 
Radius ,  cm 
Figure  4.1.   Flux  Traverses  and  Ratio  of  Lu177/ 
Lul76m  and  Pu239/Aul9?  in  i  in. 
Stainless  Steel  Rod 
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are  tabulated  in  Table  H.2   together  with  the  change 
in  effective  neutron  temperature  from  the  outer  surface 
to  the  center  line  of  the  stainless  steel  rod  using 
tabulated  g  factors.   The  ratios  show  that  the 
activity  ratio  (Lu177/1/V  absorber)  is  more  sensitive 
to  the  spectrum  shift  in  highly  absorbing  media  than 
the  (Pu239/i/v  absorber)  ratio. 

Uranium  Slab 

Lu  wires  were  placed  between  uranium  foils  as 
shown  in  Figure  4 . 2 . and  the  uranium  foils  were  positioned 
in  a  slot  in  a  graphite  block.   The  corrected  activities 
of  Lu  77  and  Lu1™*  together  with  their  ratios  are  plotted 
on  the  bottom  of  Figure  4.2  and  tabulated  in  Table  4.3. 
The  data  in  this  medium  were  interpreted  in  terms  of 
effective  neutron  temperature  using  equation  (2.8). 
Cd  ratios  were  measured  inside  and  outside  of  the  U 
slab  and  it  was  found  that  the  contribution  of  epi- 
thermal  neutrons  to  the  Lu  activities  was  negligible. 
The  effective  neutron  temperature  change  in  the  slab 
was  found  to  be  9 . 4°  C  +  3. 4°  C  using  tabulated  g 
factors. 


bJ 


.»       .6 


.i  0      .2      .4 


Figure  H.2.   Lu  Traverses  Through  U-Slab 
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TABLE  4.2 


RATIOS  OF  [Pu239/l/v  ABSORBER]  AND 
[Lu177/l/v  ABSORBER]  ACTIVITIES* 


(Pu239/l/v  absorber)       (Lu177/l/v  absorber) 
5r)+       I  C.880r)/I  (1.045r)+ 


1.000 
0.997 
0.987 
0.973 
0.948 
0.932 
0.915 
0.8977 


*     The  standard  deviations  of  these  values 
are  less  than  1%. 

t   The  change  in  neutron  temperature  using 
two  different  ratios  were: 


Radius 

I  (.980r)I  ( 

(cm) 

0.00 

1.000 

0.20 

0.999 

0.40 

0.995 

0.60 

0.989 

0.80 

0.981 

1.00 

0.971 

1.20 

0.964 

1.27 

0.957 

ATLul77   =   28.4  1  3.2°C 
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TABLE  it.  3 


Lu177  AND  Lu176m  ACTIVITIES  AND  THEIR 
RATIOS  IN  URANIUM  SLAB 


Position  (cm)     A177*      A176m*       A177/A176m* 


0.8255 

1.1300 

1.1820 

0.9560 

0.5715 

1.1134 

1.1431 

0.9740 

0.15875 

1.0126 

1.0369 

0.9766 

0 

1.0000 

1.0000 

1.0000 

-0.15875 

1.0092 

1.0267 

0.9830 

-0.5715 

1.11877 

1.1544 

0.9691 

-0.8255 

1.1360 

1.1660 

0.9743 

*   Estimated  standard  deviation  is 
less  than  1  per  cent. 


Activations  in  B2C>3  Solutions 

The  spectral  hardening  in  an  Al  can  filled  with 
aqueous  solutions  of  B203  was  measured  using  activation 
detectors.   The  detail  of  the  Al  can  is  shown  in 
Figure  4.3. 

The  Al  can,  filled  with  aqueous  B„03  solutions, 
was  placed  24  in.  from  the  east  end  of  the  thermal 
column  and  irradiated  for  ten  minutes.   The  flux  at  the 
center  of  the  can  was  approximately  4  x  109  neutrons /cm2 
sec.   Three  different  B203  solutions  were  used  as 
follows:   ordinary  tap  water,  5  grams  B^/liter,  and 
10  grams  B203/liter.   The  Maxwellian  averaged  absorption 


Figure  4.3.   Al  Can  With  the  Detector  Holder 
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cross  sections  for  these  concentrations  were  0.66  barns, 
2.609  barns  and  4.496  barns,  respectively.   The  detectors 
used  for  these  measurements  were  Dy16\  Lu175,  Lu176, 
Au197  and  Pu239. 

Normalized  activities  of  the  five  detectors  are 
listed  in  Table  4.4  and  plotted  in  Figures  4.4-4.7.   The 
ratios  of  A177/A176m  are  listed  in  Table  4.5  and  plotted 
in  Figure  4.8. 

In  order  to  compare  the  measured  activities  with 
theoretical  values,  calculations  were  made  using  the 
THERM0S  Code  (4).   The  results  of  the  calculations  gave 
the  neutron  density  and  flux  per  unit  volume  for  thirty 
energy  groups  at  different  space  positions  and  also  the 
activities  of  Dy164,  Pu239,  Lu176,  Au197  and  Eu151.   The 
calculated  activities  of  1/v  absorbers  were  found  to  be 
almost  identical  to  the  calculated  neutron  flux;  there- 
fore, experimentally  measured  activities  from  the  1/v 
absorvers  were  compared  with  the  calculated  neutron 
fluxes.   The  results  are  plotted  on  Figure  4.9  .   There 
is  some  difference  between  the  two  which  is  attributed 
to  the  inability  to  mock-up  the  axial  feed  of  the  sources 
in  the  THERMOS  code.*   The  axial  feed  of  the  sources 
makes  it  anisotropic  whereas  the  THERMOS  code  assumes 
isotropic  scattering  of  the  sources.   The  theoretical 
flux  peaks  in  the  HO  medium  a  few  mm.  inside  the  Al 

*   Kinard,  F.  E. ,  Private  Communication. 
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TABLE  4.4 
NORMALIZED  ACTIVITIES  IN  AQUEOUS  SOLUTIONS 


B203 
Cone . 

Distance 
(cm) 

Detectors 

Grams/ 
Liter 

Lu175* 

Lu176* 

Dy164* 

Au197* 

PU239* 

0 

0.6 

0.5907 

0.5955 

0.6396 

0.6386 

0.6274 

0 

1.6 

0.6969 

0.7088 



0.6733 

0.6511 

0 

2.6 

0.7814 

0.7949 

0.7679 

0.7439 

0.7374 

0 

3.6 

0.8706 

0.8731 

0.9102 

0.8655 

0.8916 

0 

4.45 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

5.0 

0.6 

0.2774 

0.2949 

0.2743 

0.2626 

0.3085 

5.0 

1.6 

0.3288 

0.3467 

0.3317 

0.3160 

0.3559 

5.0 

2.5 

0.4762 

0.5030 

0.4719 

0.4560 

0.5050 

5.0 

3.6 

0.7681 

0.7906 

0.7591 

0.7229 

0.8381 

5.0 

4.45 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

10.0 

0.6 

0.1598 

0.1788 

0.1555 

0.1604 

0.1714 

10.0 

1.6 

0.2014 

0.2227 

0.2087 

0.2061 

0.2136 

10.0 

2.6 

0.3725 

0.3986 

0.3696 

0.3439 

0.3639 

10.0 

3.6 

0.7172 

0.7402 

0.7060 

0.6607 

0.7308 

10.0 

4.45 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

Estimated  standard  deviation  less  than  1%, 
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TABLE  4.5 
A177/A176ra  RATIOS  IN  AQUEOUS  SOLUTIONS 


Distance 
Center  ( 

from 
cm) 

Light  Water 

5gm/liter* 

lOgm/liter* 

0.6 

1.0082 

1.0630 

1.0082 

1.6 

1.0169 

1.0545 

1.0169 

2.6 

1.0172 

1.0562 

1.0172 

3.6 

1.0028 

1.0293 

1.0028 

4.45 

1.0000 

1.0000 

1.0000 

*  Estimated  standard  deviation  less  than  1%. 
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can,  and  this  is  believed  to  be  due  to  the  external 
source  conditions.   The  ratios  of  A*77 '/1/v  absorber  for 
experimental  and  theoretical  calculations  are  listed 
in  Table  4.6.   The  theoretical  spectra  are  softer  than 
the  experimental  values  and  again,  this  is  due  to  the 
THERM0S  source  conditions. 


TABLE  4.6 

177 
RATIO  OF  A    /1/v  ABSORBER  IN  10  gm.  B„Cq/LITER 


Distance  (cm) Theoretical Experimental' 

0.6  1.0469  1.1191 

i-6  1.0320  1.1058 

2.6  1.0256  1.0699 

3-6  0.9991  1.0321 

4-l+5  1.0000  1.0000 


*      Standard  deviation  less  than  1  per  cent. 

In  order  to  duplicate  the  source  conditions 
imposed  by  the  THERM0S  Code,  an  attempt  was  made  to 
perform  an  experiment  by  placing  the  Al  can  filled  with 
borated  water  in  the  center  of  the  UFTR.   However,  with 
only  one  gram  of  B^/liter,  the  excess  reactivity  of 
the  reactor  was  not  enough  to  override  the  poisoning 
effect.   Therefore,  the  only  way  this  experiment  can  be 
performed  properly  is  to  provide  sufficient  reactivity 
to  override  the  B203  absorption. 
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Measurements  in  the  UFTR  Core 

In  order  to  measure  the  spacial  dependence  of 
the  effective  neutron  temperature  and  the  rvT/T0 
factor  in  the  UFTR  core,  Lu  wires  and  very  thin  Au 
foils  were  irradiated.   The  Au  foils  were  made  using 
Au  resinate  solution,*  which  was  applied  to  the  surface 
of  an  Al  foil  and  then  baked  slowly  in  an  electric  oven 
for  about  H5  minutes  until  the  temperature  reached 
900°  C.   The  baking  process  volatilized  the  organic 
material  in  the  film  of  resinate  solution  and  left  a 
thin  film  of  Au  on  the  surface  of  the  Al  foil.   The 
maximum  thickness  of  the  Au  film  was  not  more  than 
0.5  microns.   This  was  determined  by  comparing  the 
activities  of  Au  foils  of  known  and  unknown  thicknesses 
in  a  known  flux.   Since  the  self-absorption  is 
negligible  for  a  Au  foil  with  thickness  of  0 . 5  microns 
(68),  it  was  correct  to  assume  that  these  foils  were 
infinitely  dilute.   Activities  from  the  impurities  were 
found  to  be  negligible.   Accurate  weighing  of  these 
thin  Au  foils  was  practically  impossible;  therefore, 
the  same  foils  were  used  to  measure  the  bare  and  the 
Cd  covered  activities.   This  process  corrected  for  weight 


*   Supplied  by  the  Hanovia  Chemical  and  Manufacturing 
Company,  East  Newark,  New  Jersey. 
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variations.   Activity  from  the  first  irradiation  was 
subtracted  from  the  second  using  the  same  foil  for 
both  irradiations,  since  the  half-life  of  Au  is  well 
known.   The  Lu  wires  were  30  mils  in  diameter  and  con- 
tained 10  per  cent  by  weight  of  Lu  0  .   The  positions 
where  Lu  wires  and  Au  foils  were  irradiated  is  shown 
in  Figure  4.10.   One  measurement  was  made  in  the 
center  of  the  fuel  box. 

From  the  result  of  the  Au  foil  activations, 
values  of  rVT/T0  were  calculated  using  equation  (2.4). 
In  this  equation,  sQ  for  Au  is  known  (6_9)  and  is 
approximately  equal  to  17.5.   From  these  calculated 
values  of  r~V  T/TQ  the  effective  neutron  temperature  in 
the  center  of  the  UFTR  core  was  calculated  by  trial  and 
error  using  equation  (2.7).   For  these  calculations 
tabulated  values  of  g  and  s4  were  used  for  Lu176.   The 
factor  g  for  Lu   5  is  approximately  one  and  the  factor 
sQGr  was  determined  experimentally  by  the  use  of  Cd 
ratios  of  Lu175,  and  by  using  rVT/T0  obtained  from 
the  Au  foil  activation  together  with  equation  (2.5). 
The  effective  neutron  temperature  was  calculated  to  be 
76.0  t    2.9°  C  in  the  center  of  the  UFTR  when  compared 
to  the  reference  point  in  the  thermal  column. 

The  numerical  results  for  rVT/T   A177   A176m 

A177/A176m,  T  ,  R175,  R1^7  and  (s  G  )175  are  tabulated 
n   cd    cd        or 
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in  Table  4.7.   In  Figure  4.11,  A177/A176m   R197    a 

'   cd  ana 
rvT/T0  are  plotted. 

Measurements  in  the  Subcritical  Assembly 

Integral  measurements  were  made  radially  12.  S 
in.  from  the  bottom  of  the  central  fuel  tube  because  the 
slot  in  the  special  Mark  V-B  fuel  elements*  was  located 
at  this  position.   At  this  height,  the  axial  Cd  ratio 
has  reached  a  constant  value  indicating  spectral  equilib- 
rium.  An  Al  foil  holder  was  attached  to  the  tube  at 
this  height  so  that  radial  traverses  could  be  taken  to 
the  cell  boundary  and  beyond. 

Bare  and  Cd  covered  Lu  and  Au  foils  were  irradi- 
ated in  four  different  arrangements  using  two  lattice 
pitches  as  is  listed  in  Table  3.1.   Also,  THERM0S  and 
CEPTR  calculations  were  made  for  these  arrangements 
and  all  the  experimental  and  theoretical  results  are 
plotted  in  Figures  4.12  to  4.19  for  thermal  and  epi- 
thermal  activities.   All  the  experimental  results  were 
normalized  at  the  cell  boundary,  and  the  direction  of 
the  traverse  was  equidistant  between  two  fuel  tubes. 
Au  foils  0.001  in.  thick  and  0.25  in.  diameter  were  used 
together  with  foils  containing  a  10  weight  per  cent 


*   Special  unclad  Mark  V-B  and  Mark  I  fuel  elements 
with  slots  for  activation  detectors  were  borrowed 
from  the  Savannah  River  Operations  Office,  Aiken, 
South  Carolina. 
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Figure  4.14 


Au  Activity  Distribution  for  Mark 
V-B  Natural  U  With  2  2  cm  Pitch 
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Figure  4.16 


Au  Activity  Distribution  for  Mark  I  and, 
Mark  V-B  Natural  U  Fuel  With  14.7  cm  Pitch 
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Figure  4.17. 


Lu  Activity  Distribution  for  Mark  I 
and  Mark  V-B  Natural  U  Fuel  With 
14.7  cm  Pitch 
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Figure  4.18. 


Au  Activity  Distribution  for  Mark  I 
and  Mark  V-B  Natural  U  With  2  2  cm  Pitch 
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dispersion  of  Lu203  in  Al.   The  latter  were  the  0.0  30 

in.  thick  and  0.25  in.  diameter.   Activities  were 

corrected  for  the  radial  JQ(Br)  flux  distribution  in  the 

subcritical  according  to  the  relation: 

*!(R) 
*(R)   =   


where 

♦  -.(R)   =   uncorrected  flux 

*(R)   =   corrected  flux 

Epithermal  activities  were  measured  using  0.0  30  in. 

Cd  covers.   Therefore,  a  Cd  cut-off  of  0.5  ev  was  used  for 

the  THERM0S  calculations. 

In  addition,  bare  and  Cd  covered  traverses  were 
made  with  infinitely  dilute  Au  foils  to  evaluate  the 
epithermal  index  r  in  the  unit  cell  for  each  lattice 
arrangement.   Using  the  activity  ratios  of  Lu177  and 

19  8 
Au   t  effective  neutron  temperatures  at  different  cell 

positions  were  calculated  using  equation  (2.10).   For 

a  reference  the  neutron  temperature  at  the  cell  boundary 

was  determined  by  a  differential  method  in  order  to 

determine  the  effective  temperature  change  in  the  unit 

cell.   The  results  of  these  measurements  are  plotted  on 

Figures  4.20  to  4.23.   Also  shown  on  Figures  4.24  to 

4.2  7  are  the  axial  thermal  and  epithermal  fluxes  with 

the  corresponding  Cd  ratios. 
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Figure  4.24. 
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CHAPTER  V 

RESULTS  OF  DIFFERENTIAL  SPECTRUM  MEASUREMENTS 

Introduction 

Differential  spectrum  measurements  were  made 
using  the  crystal  diffraction  spectrometer  described 
previously.   A  beam  of  neutrons  was  extracted  from  the 
subcritical  assembly  for  this  purpose.   This  choice  of 
instrument  was  made  because  the  crystal  spectrometer 
is  simpler,  cheaper  and  less  space-consuming  than  a 
chopper  (70).   The  details  of  the  crystal  spectrometer 
and  the  description  of  the  subcritical  assembly  are 
given  in  Chapter  III. 

Analysis  of  the  crystal  spectrometer  data  is 
complicated  due  to  the  following  effects:   (1)  secondary- 
and  multiple  scattering  in  the  walls  of  the  collimator 
(2)   the  Renninger  effect  (71)*   (3)  imperfections  of 
the  crystal,  i.e.,  mosaic  structure  (H)      variations  in 
detector  efficiency,  and   (5)  higher  order  contamination. 
In  order  to  simplify  these  corrections,  the  instrument 


The  Renninger  effect  is  explained  in  detail  in 
Chapter  III. 
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was  calibrated  experimentally  by  using  a  beam  of  neutrons 
with  a  known  energy  spectrum.   More  detailed  discussion 
of  this  method  is  given  in  the  later  sections  of  this 
chapter.   For  calibration  purposes,  a  well-moderated 
neutron  beam  was  extracted  from  the  subcritical  tank 
filled  with  D20. 

Following  this  calibration  experiment,  beams 
from  the  D20-moderated  subcritical  were  extracted  from 
the  center  of  the  fuel  element,  from  the  fuel-D20 
interface  and  from  the  moderator  surrounding  the  fuel 
elements.   A  similar  work  was  performed  by  V.  I.  Mostovoi, 
et  al.  (7^)  in  uranium-water  and  uranium  monoisopropyl- 
diphenyl  lattices  using  the  time-of- flight  technique. 

For  theoretical  comparisons,  the  THERM0S  code 
developed  by  Honeck  <»0  was  used.   THERM0S  is  a  multi- 
group  code  which  solves  the  integral  form  of  the 
transport  equation;  it  is  described  in  Appendix  B  . 

Alignment  and  Calibration 
Alignment 

Several  different  methods  of  aligning  were  tried. 
First,  the  bore-sight  method  of  aligning  the  detector 
and  the  first  collimator  was  tried.   Then,  an  optical 
method  was  attempted.   Finally,  it  was  decided  that  the 
easiest  way  to  align  the  spectrometer  was  to  set  the 
detector  shield  in  the  vertical  position  and  then  to 
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insert  an  Al  rod  1  3/8  in.  diameter  through  the  detector 
shield  opening  with  guides  attached  to  the  bottom  and  top 
of  the  shield.   Then,  a  guide  to  the  top  of  the  collimator 
was  attached  so  that  when  the  detector  hole  and  the 
collimator  were  aligned  very  accurately,  the  Al  rod 
would  fit  into  the  guide  on  top  of  the  collimator.   The 
Al  rod  was  within  ±0.001  in.  from  being  straight.   In 
fact,  a  deviation  of  ±0.001  in.  in  the  Al  rod  will  cause 
an  uncertainty  in  6  by  ±  one  minute.   The  final  method  of 
aligning  is  shown  in  Figure  5.1. 

Spectrum  From  the  UFTR  Core 

In  order  to  measure  the  differential  neutron 
spectrum  in  the  center  of  the  UFTR  core  the  central 
vertical  plugs  above  the  reactor  were  replaced  with 
a  collimator  plug.   The  diffraction  spectrometer  was 
then  aligned  with  the  collimator.   Rocking  curves 
obtained  with  NaCl  (200)  and  LiF  (111)  crystals  are 
plotted  on  Figures  5.2  and  5.3.   Half  maximum  angles, 
A 6,  of  the  NaCl  and  LiF  were  1.7°  and  0.5  2  5° 
respectively. 

Using  the  LiF  crystal,  the  neutron  spectrum 
coming  out  of  the  central  vertical  access  was  measured 
from  3°  <  e  i  27°.   Data  obtained  from  this  measurement, 
after  correcting  for  the  background,  were  plotted  as  E 
versus  count  rate  as  shown  in  Figure  5.H.   It  can  be 
seen  from  the  figure  that  the  results  fitted  a 
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Maxwellian  distribution  very  well.   Therefore,  the  neutron 
temperature  was  calculated  for  the  Maxwellian  distribution. 

According  to  the  kinetic  theory  of  gasses  (73), 
the  average  kinetic  energy  of  thermal  neutrons  is  given 
by: 

i  mv2   =  3   k  Tn 

2         2 

(5.1) 

The  quantity  v2  is  not  measurable  but  it  is  related  to 
vQ,  the  most  probable  velocity,  which  is  measurable  as 
shown  below  ( 74) . 


v* 
/v2dn      aAvq3  Jv"e~  v?  dv 


(5.2) 


7 

/dn  H  -   2  vo' 


The  most  probable  velocity,  vQ,  of  the  Maxwell  dis- 
tribution is  obtained  by  differentiating  the  flux 
distribution  (equation  (5.3))  and  setting  it  equal 
to  zero. 


d(nv)       4n         _  JL, 
dv    "   -^v3   v3  e   v02  ^^ 


or 


nu   =   1.224  v 


(5.4) 
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where  vm  corresponds  to  the  velocity  of  the  most 
probable  flux  or  root  mean  square  velocity.   In  a 
simpler  manner,  the  neutron  temperature  of  the  core  is 
related  to  the  most  probable  velocity  as  follows: 


irav„2   =   k  T 
o     o         n 


(5.5) 


The  neutron  temperature  for  Figure  4.9  was 
calculated  as  outlined  above  and  was  found  to  be 
115.4  t   6.0°  C  with  the  most  probable  velocity  of 
2531  t    314. 5  m/sec. 

Spectra  From  D20 

In  an  infinite  moderator  with  vanishingly  small 
absorption,  the  neutron  spectrum  takes  the  form  of  a 
Maxwellian  distribution  at  the  temperature  of  the  media. 
If  a  poison  with  1/v  absorption  cross  section  is  added 
to  a  moderator,  it  can  be  shown  (75_)  that  if  the  Max- 
wellian spectrum  is  to  be  retained,  the  source  of  the 
neutrons  must  itself  be  a  Maxwellian  spectrum.   In 
this  case,  the  source  and  loss  rate  from  each  energy 
interval  must  be  the  same.   In  an  actual  reactor,  the 
spectrum  is  distorted  from  Maxwellian,  and  this  is 
due  to  the  slowing  down  distribution  of  the  source 
from  fission,  which  does  not  have  the  Maxwell  distribution. 
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A  neutron  beam  has  been  extracted  from  scatterers 
such  as  H20,  D20  and  C  (76^  TJ_)    and  the  neutron 
spectrum  analyzed  with  the  time-of-flight  method  using 
a  chopper  and  neutron  pulser.   The  measurements  confirmed 
that  the  neutron  energies  fit  a  Maxwellian  distribution 
very  v/ell  from  the^e  non-abqorbina  media. 

In  order  to  extract  a  beam  of  neutrons  from 
non-absorbing  media,  a  cylindrical  tank  2U  in.  diameter 
and  60  in.  high  was  placed  on  top  of  the  graphite 
pedestal  of  the  UFTR  and  was  filled  with  D20.   a 
more  detailed  description  of  the  set-up  is  discussed 
in  Chapter  III. 

The  purpose  of  extracting  a  beam  of  neutrons 
out  of  the  D2O  medium  with  variable  temperature  is  to 
calibrate  the  overall  efficiency  of  the  neutron  diffrac- 
tion spectrometer  using  neutron  spectra  with  fairly  well- 
known  distributions. 

The  count  rate  of  a  BF   counter  at  different 
Bragg  angles  6  for  a  single  crystal  can  be  expressed 
as: 


n=k 
Ci<6)   =    I      TifnFi.n(6) 


n=1  (5.6) 
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where 

T^  =  transmission  of  the  nth  order 
diffracted  neutrons  through  a 
filter  (e-EtX) 


i  =  measurement  index 

n  =  order  of  diffraction 

C^(8)  ■  count  rate  at  the  detector 

F.  „(e)  =  *(n2E)e(n2E)S(n2E)  w(n2E)  r(n2E)  D(n2E) 
ii"  n      n      n     n 

<|>(n2E)  =  neutron  flux  distribution  in  the  media 


e(nzE)   =   BF3  detector  efficiency 

2 
S(n  E)    =   spectrometer  resolution  for  order  n 

2 
W(n  E)n   =   fraction  of  reflected  beam  transmitted 

by  the  collimating  system  for  order  n 


r(n  E)n   =   reflectivity  of  the  crystal  plane  used 
for  order  n 


D(n  E)n   =   correction  factor  due  to  the  multiple 
diffraction  for  order  n 


The  evaluation  of  F.  _(8)  in  equation  (5.6)  was 
x  ,n 

carried  out  for  up  to  two  sets  of  T^  n's,  by  solving 

simultaneous  equations.   Since  the  $(n2E)  is  known  for 

large  D20  media,  the  combined  efficiency   for  different 

energies  for  first  order  can  be  evaluated.   Instead  of 

evaluating  the  combined  efficiency,  ratios  of  F.   (e) 

l  ,n 
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of  unknown  and  known  spectra  were  taken  and,  when  multiplied 
by  the  known  spectrum,  gave  the  unknown  energy  spectrum: 


♦<n2E)x   =   ,(n2E) 

CFi  n^e>^Q  S 

x»n     s  (5.7) 


The  neutron  spectrum  was  extracted  7  in.  from  the 
bottom  of  the  tank  with  D20  temperatures  of  25°  and  75°C. 
The  axial  Cd  ratio  for  this  system  is  plotted  on  Figure 
4.27.   By  inspection  of  the  data  it  was  found  that  the 
neutron  temperature  was  10° C  higher  than  the  moderator 
temperature. 

The  ratio  of  the  two  sets  of  measurements  are 
plotted  on  Figure  5.5  with  the  theoretically  calculated 
Maxwellian  distribution  ratios  at  two  different  temperatur 
This  was  done  to  check  if  the  spectrometer  could  see  the 
change  in  neutron  spectrum  due  to  the  change  in  medium 
temperature. 

Spectra  from  Subcritical  Assembly 

For  the  differential  spectrum  measurements,  the 
first  three  lattice  arrangements  listed  in  Table  3.1  were 
used.   The  first  collimator  described  earlier  was  inserted 
vertically  into  the  subcritical  assembly  first  at  the 
center  and  then  at  the  boundary  of  the  central  unit  cell. 
The  LiF  (111)  crystal  was  used  for  these  measurements. 
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Figure  5.5. 


Ratio  of  Maxwellian  Distribution  and 
the  Experimental  Spectrum  at  Two 
Different  Temperatures 
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The  experimental  neutron  density  distributions  were 
calculated  at  two  locations  by  the  use  of  equation  (5.7) 
for  each  arrangement  and  are  plotted  on  Figures  5.6  -  5.8. 
In  Figure  5.7  the  lower  experimental  spectrum  is  much  softer 
than  the  THERM0S  spectrum.   However,  the  THERM0S  spectrum 
is  at  the  center  of  the  Mark  I  fuel  element,  while  the 
extracted  beam  of  neutrons  had  a  diameter  of  approximately 
5  cm.   Therefore,  the  experimental  results  were  plotted  again 
in  Figure  5.9  with  THERM0S  spectrum  0 . 2  cm  outside  of  the 
Mark  I  fuel  element. 
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Mark  V-B  Natural  U  Fuel  With  2  2  cm  Pitch 
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Figure  5.7.   Experimental  and  THERM0S  Spectra  Using 
Mark  I  and  Mark  V-B  Natural  U  Fuel 
With  22  cm  Pitch 
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Figure  5.8.   Experimental  and  THERM0S  Spectra  Using 

Mark  V-B  Natural  U  Fuel  With  14.7  cm  Pitch 
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Figure  5.9.   Experimental  and  Corrected  THERM0S 
Spectra  Using  Mark  I  and  Mark  V-B 
Natural  U  Fuel  With  22  cm  Pitch 


CHAPTER  VI 

RESULTS,  CONCLUSIONS  AND  RECOMMENDATIONS 

Comparisons  of  Experimental  and  Theoretical 
Spectra  in  the  UFTR  Thermal  Column  and  Core 

Comparisons  of  the  theoretical  activities  with 
the  experimental  results  in  the  thermal  column  using  an 
Al  can  filled  with  borated  water  are  plotted  on  Figure 
4.9.   The  differences  in  the  two  fluxes  were  as  much  as 
2  5  per  cent  and  this  is  due  to  the  THERM0S  source  con- 
ditions.  It  would  be  interesting  to  perform  this  experi- 
ment with  isotropic  source  conditions  in  a  reactor  which 
has  enough  excess  reactivity  to  override  the  poisoning 
effect  of  the  sample. 

The  ratio  of  the  activities  of  a  resonance  and  a 
1/v  absorber  for  borated  water  was  found  to  be  somewhat 
higher  than  the  theoretical  calculations  as  can  be  seen 
in  Table  4.6.   The  discrepancy  is  again  due  to  the  THERM0S 
source  conditions. 

No  theoretical  calculations  were  made  for  the 
natural  U  slab  and  the  stainless  steel  rod.   Again,  anisotropy 
of  the  source  conditions  would  apply  here  and  there  would 
probably  be  large  discrepancies  between  THERM0S  calculations 
and  experimental  results  as  with  the  borated  water  data. 
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From  the  experimental  results  a  good  comparison  was 
made  between  the  integral  and  the  differential  measurements 
in  the  center  of  the  UFTR  core.   The  effective  neutron 
temperature  determined  by  the  differential  measurement 
was  115±6°C  compared  to  76±2.9°C  using  integral  measurements 
with  Westcott's  formulations  (3H_).   Thus,  the  difference 
between  the  two  methods  of  measurements  was  39° C.   The 
effective  temperature  of  the  integral  measurement  in  the 
UFTR  core  is  relative  to  the  effective  temperature  of  the 
thermal  column  which  was  assumed  to  be  equal  to  the 
temperature  of  the  media.   Johansson  (]J>)    nas  measured  the 
effective  neutron  temperature  in  graphite  media  using  the 
time-of-f light  method  and  found  it  to  be  29±10°C  higher 
than  the  medium  temperature.   If  this  result  applies  to 
the  UFTR  thermal  column,  then  the  comparison  of  effective 
temperature  using  the  two  methods  will  be  within  the  experi- 
mental error. 

The  experimental  effective  neutron  temperature  in 
the  UFTR  was  compared  with  theoretical  calculations  made  by 
C.  A.  Thompson  (7jO  who  has  made  a  multigroup-multiregion 
calculation  for  one  dimension  going  from  South  to  North 
in  the  reactor  core.   However,  his  calculations  showed  a 
much  lower  neutron  temperature  in  the  UFTR  core.   It  is 
also  hard  to  make  a  detailed  comparison  with  Thompson's 
calculations  because  his  energy  groups  in  the  thermal 
region  are  wide. 
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Comparisons  of  Experimental  and  Theoretical 
Spectra  in  the  Subcritical  Assembly 


Intracell  activation  measurements  were  made  using 
resonance  and  1/v  absorbers  in  each  of  the  lattice  arrange- 
ments in  Table  3.1,  and  comparisons  were  made  with  THERM0S 
calculations  for  each  detector.   In  general,  the  results 
are  in  good  agreement  except  for  a  few  points  like  the 
one  in  Figure  4.14.   These  few  points  are  very  probably 
due  to  inaccurate  positioning  of  the  activation  detectors 
in  the  horizontal  plane.   This  conclusion  was  verified  by 
repeating  some  of  the  measurements.   Also,  CEPTR  calculations 
for  the  flux  distribution  were  made,  but  were  found  to  be 
almost  identical  to  the  THERM0S  flux. 

For  the  differential  measurements ,  experimental 
calibration  of  the  equipment  proved  to  be  successful  and 
gave  excellent  results,  as  can  be  seen  in  Figures  5.6  - 
5.9.   The  temperature  calibration  curve  was  also  satisfactory 
which  can  be  seen  in  Figure  5.5. 

The  comparisons  of  the  change  in  effective  neutron 
temperature  from  cell  boundary  to  cell  center  for  three 
different  methods  (THERM0S  calculation,  integral  and 
differential)  are  listed  in  Table  6.1. 

The  temperature  change  in  a  unit  cell  for  the 
differential  measurements  was  estimated  by  fitting  the 
experimental  3pectrum  to  a  Maxwellian  distribution  using  a 
least  square  fit  with  the  IBM  704  computer  ( 81) . 
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Ill 

It  would  be  of  interest  to  extend  the  data  to 
other  fuel  element  configurations  (e.g.,  large  clusters 
of  U02  rods)  in  order  to  test  the  validity  of  the  THERM0S 
calculations. 


APPENDIX  A 
METHOD  OF  CALCULATING  A177/A176m 

Natural  Lu  contains  2.6  per  cent  Lu176  and 
97.4  per  cent  La*7*.   Activities  of  irradiated  natural 

Lu  comes  from  Lu175m  and  t,,177  .   ..   ...   .  . 

ana  lu    Xn   the  fxrst  two  days. 

After  four  days  the  activity  of  Lu176m  decays  out 

completely,  and  only  the  activity  of  Lu177  is  left. 

Normally,  each  foil  or  wire  was  counted  ten  times  to 

get  better  statistics.   The  first  five  counts  were  made 

for  the  short  and  long  half-life  activities,  and  then 

another  five  counts  were  taken  for  the  long  life 

activity.   All  the  activities  were  corrected  back  to 

some  standard  time,  t=0,  and  activities  normalized  to 

standard  Cs137  source. 

The  average  activity  of  Lu177  at  t=0  is  given  as: 

A  t. 

5   [A*77  (tjL)  -  Background]  e  2  1 

A*77  (t  =  0)   =    I      — +  5177 


i  =  l 


x 


where 


decay  constant  of  Lu177 

time  when  activity  of  Lu176m  and 

17  7 
Lu    was  measured 
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x   «   index  for  the  position  of  the 
detector 

5^    =   standard  deviation  of  Lu177 
activity 

Now  the  average  activity  of  Lu176m  at  t=0  is  given  as 


A176m 

Aj     (t=0) 


5  .  . 

.  rA176m  +  177,      177         2  j  Xltj 

I    LAX         (t.)-A   (t=0)  e      -  background]  e 
j  =  l  J    x 


+   ,176m 
"   6x 


\1      -      decay  constant  of  Lu176m 

t.   =   time  when  the  activities  of 
]      Lu?76m  and  Lu177  were  measured 

6^76m  =   standard  deviation  of  the  Lu176m 
activity 


All  the  activities  were  normalized  to  standard  activities, 

A177(t=0) 
Y177(t  =  0)   =  — ,   A177, 


*?V"        ~ -   Ax"< 

A^77(t=0) 
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A^75m(tS0) 
^76m(t  =  0)   c   s   Ai76m(t  =  0) 


A^76m(t=0) 


x 


Since  A^   =1   and  Ai/bm  =  1.   The  normalized  activity 
ratios  become 


Al?7(t=0) 
r   =  ±   .  R 


A176m(t=0) 


*   6 


x 


APPENDIX  B 
THE  THERM0S  CODE 

The  THERM0S  Code  developed  by  H.  C.  Honeck  of 
Brookhaven  National  Laboratory  calculates  the  scaler 
neutron  spectrum  as  a  function  of  energy  and  position 
in  a  lattice  cell.   The  code  solves  numerically  the 
integral  transport  equation  with  isotropic  scattering  (4). 
One  dimensional  slab  or  cylindrical  geometry  can  be  used. 
Comparisons  of  the  code  calculations  with  integral  measurements 
have  been  published  (79,  80). 

THERM0S  code  uses  30  energy  groups  as  given  in 
Table  B.l;  it  uses  20  space  points  and  mixtures  of 
ten  nuclides.   The  code  is  used  for  calculations  of  the 
reaction  rate  of  11*33,  u235  §  Pu239(  Dy16^  Eu151f  and 
Lu176.   Output  of  the  THERM0S  code  is  also  used  for 
calculating  n,f  (thermal  utilization),  £  ,  fQ  (fraction 
of  thermal  neutrons  captured  in  U238)  and  D. 
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TABLE  B.l 
ENERGY  GROUPS  FOR  THERM0S  CALCULATIONS 


Group  No. 

Energy  (ev) 

Group  No. 

Energy  (ev) 

1 

0.00025 

16 

0.06517 

2 

0.00101 

17 

0.07485 

3 

0.00228 

18 

0.08612 

4 

0.00405 

19 

0.09919 

5 

0.00532 

20 

0.11398 

6 

0.00911 

21 

0.13123 

7 

0.01240 

22 

0.1525 

8 

0.01619 

23 

0.1790 

9 

0.02049 

24 

0.2124 

10 

0.02530 

25 

0.2546 

11 

0.03061 

26 

0.3081 

12 

0.03643 

27 

0.3766 

13 

0.04276 

28 

0.4618 

14 

0.04959 

23 

0.5402 

15 

0.05692 

30 

0.7067 

APPENDIX  C 

THE  RECIPROCAL  LATTICE 

Let  a-,  ,  a2  and  a3  represent  vectors  of  a  unit 
cell  containing  no  lattice  points  other  than  those 
lying  at  its  corners.   Then,  let  blt  fc>2  and  t>3  represent 
the  vectors  of  a  reciprocal  lattice.   The  reciprocal  axis 
b-,  is  perpendicular  to  the  axis  a2  and  a3.   Also,  b2 
is  perpendicular  to  a-j^  and  a3  and  t>3  is  perpendicular 
to  a-,  and  a2 .   These  relations  are  summed  up  in  the 
vector  equations: 

tH  .  a0  =  bn  •  a^  =  b0  .  a-,  =  b0 


(C.l) 


The  magnitude  of  reciprocal  vectors  are  fixed 
by  relations: 


bl  "  d2  '  b2  S  a3  *  b3  =  1  (C'2) 


The  length  of  the  reciprocal  lattice  is  given 
by  relation: 


b1    =  3(42«  a3)  (C.3) 
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where  B  is  a  constant  evaluated  as  follows :   Let  V  be 
the  volume  of  the  unit  cell  of  the  lattice  defined  by 
vectors  a  ,  a  ,  a  ,  then: 


<a2  x  a3)  =  a2 


Ca3  x  ax) 


(ax  x  a2) 


(C.4) 


If  y  is  the  angle  between  £  and  a2,  then  in  Figure  C.l 

al  x  a2 '  r  'a'  'b'  sin  y   is  the  area  of  the  face 
OA,A4A2  so  that  a3  •  (a,  x  a2)  =  a3  cos  6  (area  of  the 

base  of  the  cell)  and  is  thus  equal  to  the  volume  V. 
From  equations  (C.3)  and  (C.U  )t   the  following  relation- 
ship is  obtained. 


V   =   (ax  •  b1)/6   =   l/l 


(C.5) 


b,  =  (a2  x  a3>,  b2  =  (a,  x  a-,),  b,  =  (a,  x  a2 ) 


(C.6) 


Figure  C.l.   Unit  Cell  of  a  Lattice 
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The  vector  r(u,v,w)  to  any  lattice  point  is  given 
by  vector  equation: 


r(u,v,w)   =   ua,  +  va2  +  wa- 


(C.7) 


where  u,  v,  and  w  are  any  integers,  positive,  negative 
or  zero.   The  vector  r'(h,k,i)  to  any  point  in  reciprocal 
lattice  is  given  by: 

r'(h,k,i)  =  hb^  +  kb2  +  £b3  (C.8) 

where  h,  k,  I   have  the  same  properties  as  u,  v,  and  w. 
Two  properties  of  the  lattice  that  make  it  of  great  value 
in  treating  problems  of  diffraction  by  space  lattices  are: 

(1)  The  vector  r'(h,k,t)  to  the  point  (h,k,X.) 
of  the  reciprocal  lattice  is  normal  to  the 
planes  (h,k,£)  of  the  crystal  lattice. 

(2)  The  magnitude  of  the  vector  r'(h,k,l)  is 
the  reciprocal  of  the  spacing  of  the  planes 
(h,k,£)  of  the  crystal  lattice. 


APPENDIX  D 


BUCKSH0T  CALCULATIONS 


CEPTR  fluxes  are  used  as  input  to  the  BUCKSH0T 
calculations.   It  calculates  the  lattice  parameters  of 
natural  uranium  D  0- moderated  cylindrical  heterogeneous 
reactors.   Some  of  the  parameters  obtained  for  two 
different  pitches  using  Mark  V-B  elements  are  listed 
below. 


•itch 

14.7  cm 

2  2.0  cm 

L2 

41.7  cm2 

115.6  cm2 

M2  2 

180.6  cm2 

2  4  3.1  cm2 

P 

0.765 

0.887 

k 

00 

1.024 

1.181 

R  2 
Bm 

1.308xl0-4  cm'2 

7.030xl0"4  cm" 

n 

1.327 

1.327 

T 

13  8.9  cm2 8 

12  7.5  cm2 

f 

0.981 

0.975 

e 

1.029 

1.029 
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APPENDIX  E 


PROPOSAL  TO  USE  UFTR  AS  A  SOURCE  FOR 
THE  SUBCRITICAL  ASSEMBLY 


PROPOSAL  August  3,  1964 

TO:      UFTR  Subcommittee  of  the  ENOF  Committee 

FROM:    S.  Salah 

SUBJECT:  Use  of  UFTR  as  a  Source  for  a  2  4"  D  x  60" 
Subcritical  Reactor 

Description  of  the  Set-up 

A  series  of  experiments  is  planned  with  a 
subcritical  reactor  sitting  on  top  of  the  UFTR  as 
shown  in  Figure  1.   A  new  set  of  shielding  blocks 
with  heavy  aggregate  will  be  constructed  which  will 
be  used  as  permanent  shielding  blocks  for  the  top 
of  the  reactor. 

When  the  subcritical  is  removed  from  the 
shield  structure,  it  will  be  replaced  with  permanent 
plugs  with  vertical  access  holes.   The  shielding 
plugs  will  be  made  from  heavy  aggregate  concrete  and 
will  be  lifted  out  with  the  crane  whenever  necessary. 

Object  of  the  Experiment 

The  purpose  of  these  experiments  is  to  measure 
the  differential  and  integral  spectrum  emerging  from 
the  subcritical  reactor.   These  measurements  will 
comprise  three  parts: 

(1)   Using  submerged  heaters  in  the 
subcritical  filled  with  H-O  and 
then  with  D20,  a  calibration 
experiment  will  be  done  by  heat- 
ing the  media.   Following  this, 
the  same  experiments  with  borated 
H20  will  be  performed. 
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(2)   The  subcritical  with  a  lattice 
spacing  of  5.8  in.  and  8.64  in. 
using  Mark  V-B  type  fuel  elements 
will  be  used.   A  beam  will  be 
extracted  from  the  fuel  element 
and  from  the  D20  moderator  and 
hardening  of  the  neutron  spectrum 
in  the  fuel  will  be  investigated 
using  a  crystal  diffraction  spec- 
trometer and  activation  foils. 
The  moderator-to-fuel  ratios  of 
the  two  lattice  spacings  are 
12.1  and  24.8,  respectively. 


Safety  Measures 

The  coupling  between  the  subcritical  and  the 
critical  reactor  will  be  measured  by  noting  the  critical 
rod  settings  before  and  after  the  installation  of  the 
subcritical  reactor.   For  the  first  reaotor  run  following 
the  installation,  the  usual  start-up  procedures  for  a 
new  fuel  loading  will  be  followed. 

At  very  low  power  level  the  flux  level  in  the 
subcritical  will  be  measured  and  compared  with  some 
reference  point  in  the  reactor.   The  area  above  the  top 
of  the  subcritical  will  be  monitored,  and  if  the  radiation 
level  is  not  above  the  permissible  limit,  it  is  proposed 
to  take  the  power  level  of  the  UFTR  to  the  maximum  level 
for  the  purpose  of  extracting  a  beam.   A  shielding  block 
will  be  placed  on  top  of  the  subcritical  for  100   kw 
operation  if  it  is  found  to  be  necessary. 

Flux  levels  and  build  up  of  activities  in  the 
central  fuel  element  is  estimated  in  Appendix  A.   After 
the  experiment,  the  lattice  of  the  subcritical  tank  will 
be  lifted  out  of  the  reactor  shield  and  it  will  be  stored 
in  the  water  shield  tank  of  the  UFTR  until  the  radio- 
acitvity  of  the  fuel  elements  goes  below  the  permissible 
limit. 

A  request  for  performing  this  experiment  has  been 
approved  by  the  Atomic  Energy  Commission. 
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Appendix  A 


I.   Calculation  of  the  flux  in  the  center  line  of 
the  subcritical  assembly. 


The  general  time  independent  one-velocity  age 
diffusion  equation  for  multiplying  media  for  cylindrical 
coordinates  is  expressed  as 


13     3  32 

[r  —  ♦(r,  Z)]  +  "-S-  . 

r  3r    3r  3Z2  m  (1) 


[r  7-  ♦(r,  Z)]  ♦  7^2"  *(r»  Z)  +  Bm  *(r»  z)  =  ° 


where 

„2    In  k« 


1/  +  T 


The  solution  for  the  flux  in  the  center  line  of  the 
cylinder  is  given  by 

♦(Z)  =  C  e"aZ  [1  -  e'2a(H  "  Z)]  (2) 


where 


Z  =  0  at  the  reactor  pedestal  subcritical 
boundary  line 


Pf 


2     2 


H  =  extrapolated  height  of  the  subcritical 
assembly. 


Four  factors  of  k„  were  calculated  by  using  the  method 
outlined  in  Murray,  "Introduction  to  Nuclear  Engineering,' 


124 
Chapter  9  and  the  following  values  were  obtained: 

n  =  1.32 
f  =  0.994 
P  =  0.876 
c  =  1.003 
k»  =  1.152 


L2  =  L2  (1-f)  =  1.346  x  104  (.006)  =  80.74  cm2 
o 


t  =»  12  5  cm' 


Bm2  =  6.9  x  10-4  cm-2 


which  gives   a  =  0.082  cm 

Using  Eq.(2)  the  flux  can  be  calculated, 


II.    Dose  from  the  irradiated  fuel  elements  in  the 
center  of  the  subcritical  assembly. 

The  total  number  of  fission/sec  in  the  fuel  element 

can  be  calculated  by  evaluating  the  following  integral: 


Z  =  100 

-aZ     p 

C  e     A  If  dZ 


L 


(3) 


X 


125 

where 

[1  -  e-(H  "  Z)] 

term  can  be  neglected  in  the  above  equation. 
C  «  <t>0  evaluated  at  Z  =  0 . 

A  =  surface  of  the  annular  part  of  the  fuel 
element . 

At  reactor  power  level  of  100  kw,  4>Q  is  approximately 
10l0  n/cm2-sec,  which  was  extrapolated  from  "Guide  to 
Irradiations  in  the  University  of  Florida  Training 
Reactor"  by  A.  R.  Boynton. 

The  integral  was  evaluated  in  this  power  level 
and  the  number  of  fissions  per  sec  was  equal  to  3.32 
10H.   This  number  corresponds  to  10.1  watt-sec. 

Irradiation  for  10  hours  at  this  power  level 
gives  a  decay  rate  of  the  fuel  element  of  3.03  x  10-L-L 
Mev/sec  one  hour  after  reactor  shutdown,  which  was 
obtained  from  Figure  3-12  of  "Fundamental  Aspects  of 
Reactor  Shielding"  by  Goldstein.   Fifteen  hours  after 
shutdown  of  the  reactor  the  activity  goes  down  to 
4.55  x  10i0  Mev/sec. 

Since  the  activity  of  the  fuel  element  will  be 
hottest  on  the  lower  20  cm  of  the  fuel  element,  the 
dose  rate  of  60  cm  from  the  lower  part  of  the  fuel 
element  was  calculated  for  100  kw-10  hour  operation 
using  formula  (4)  out  of  "Reactor  Shielding  Design 
Manual"  by  Rockwell: 

S  R  2 
v  o 

D(Mev/cra2-sec)  = F(6,  b2) 

2(a  +  Z) 


and  found  to  be  2000  mr/hr  at  15  hours  after  reactor 
shutdown,  or  100  mr/hr  after  26  days. 


(4) 


126 
III.    Total  Pu239  build-up. 

120  cm 


,239    r      .  «  _"oZ  _  2  3 : 

0 


:  f      <t'oA  e     ^J" 


1010  x  13.61 

~~      x  .1286  x  3.6  x  10H 
.082 


=  7.68  x  1015  Pu239 


atoms 


or,  weight  Pu239  =  39  y  grams 


IV.    Total  fission  product  build-up, 


120  cm        qZ 
FP  =  f  ♦Q  A  e"    z235   dZ  (total  time) 

Jo 


f 


where 

fission/sec  =   3.32  x  1010 

fissions  in  10  hrs .   =   3.32  x  10^  x  3.6  x 


104  =  7.68  x  10 
Pu239  atoms 


IS 


number  fissions      =   1.55  x  10 
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